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SECTION  I 
INTRODUCTION 

For  some  time  now,  the  AFFOL  has  been  investigating  the  use  of 
rotary-vaned  positive-displacement  air-cycle  machines  for  aircraft 
environmental  control  systems.  The  selection  of  a rotary-vaned  unit 
is  based  on  the  assumption  that  a positive  displacement  unit  of  this 
type  is  capable  of  operating  more  efficiently  than  conventional  turbo- 
machinery  at  the  low  flow  rates  required  to  cool  distributed  small  heat 
loads. 

Preliminary  rotary  cooler  studies  were  undertaken  with  two  major 
objectives  in  mind:  (1)  to  identify  the  most  promising  combination  of 
rotary  machine  and  equipment  arrangement  for  meeting  the  cooler  require- 
ments, and  (2)  to  identify  any  significant  problem  areas  and  possible 
solution  approaches. 

Two  types  of  rotary  positive-displacement  machinery  were  considered 
during  previous  investigations;  they  were  the  single-rotor  sliding-vane 
type  of  machine,  and  the  two-rotor  fixed-vane  machine.  Although  there  are 
a number  of  different  machinery  designs  which  use  the  basic  sliding-vane 
principle,  only  one  design  configuration  (herein  referred  to  as  the 
ROVAC  machine)  was  considered  for  the  rotary  cooler  application.  Likewise, 
only  one  configuration  of  the  two-rotor  fixed-vane  machine  was  considered. 
These  two  machines  had  considerably  different  mechanical  features  and 
thus  provided  a reasonable  basis  for  initial  assessment  of  machinery 
viability  for  the  rotary  cooler  application. 

Ultimately,  a sliding-vane  or  ROVAC  machine  was  selected.  Schematic 
drawings  of  the  ROVAC  sliding-vane  machine  are  shown  in  Figures  1 and  2. 

The  machine  consists  of  a circular  slotted  rotor,  concentrically 
positioned  in  the  elliptical  bore  of  a stationary  housing.  Sliding 
vanes  in  the  rotor  slots  divide  the  space  between  rotor  OD  and  housing 
bore  into  a number  of  small  spaces  (cylinders)  whose  volumes  vary  as 
the  rotor  turns.  On  one  side  of  the  bore's  minor  axis  these  variable 
volumes  are  used  to  accomplish  positive  displacement  compression;  on 
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the  other  side  of  the  minor  axis  the  variable  volumes  perform  positive 
displacement  expansion.  As  depicted  in  Figure  2,  cams  and  cam  rollers 
at  the  two  ends  of  the  rotor  can  be  used  to  actuate  the  sliding  vanes 
within  the  rotor  slots.  Or,  vane  actuation  may  be  achieved  via  contact 
forces  between  the  vane  tips  and  the  cylinder  wall,  as  is  usually  the 
case  in  commercial  sliding-vane  machines.  The  cam  method  of  actuation 
permits  reduction  or  elimination  of  vane-tip  contact  forces,  thus  reducing 
or  eliminating  vane-tip  friction  losses  and  wear.  The  vane-tip  method 
of  actuation  is  mechanically  simpler  (and  lower  in  manufacturing  cost) 
and  will  usually  be  selected  for  applications  where  vane-tip  friction  and 
wear  do  not  present  severe  problems. 

The  ultimate  method  selected  for  engineering  development  was  the 
cam  actuation  method  because  of  its  potentially  higher  efficiencies. 
However,  this  method  has  inherent  problems  in  mechanical  complexity. 

Within  the  engineering  development  the  true  severity  of  the  mechanical 
problems  associated  with  cam  actuation  became  apparent.  These  problems 
are  especially  true  at  the  even  higher  design  speeds,  where  payoffs  in 
increased  efficiencies  and  reduced  weight  become  highly  desirable. 
Basically,  two  distinct  problems  presented  themselves  and  consisted  of: 

(1)  bearing  load/speed  limitations  and  (2)  excessive  vane-tip  deflection 
values.  Figure  3 illustrates  both  of  these  problems  on  a typical  vane 
for  use  in  a cam-actuated  compressor. 

The  normal  solution  to  either  one  of  the  two  problems  inherently 
leads  to  compounding  the  other.  For  example,  vane-tip  deflection  can  be 
reduced  via  heavier  vanes;  however,  increased  bearing  loads  will  occur 
due  to  higher  centrifugal  forces  from  the  increased  mass.  Conversely, 
reduced  vane  masses  generally  reduce  bearing  loads  but  lead  to  increased 
tip  deflections. 

One  potential  solution  to  the  load/deflection  problem  is  the  use 
of  advanced  low-density  graphite  materials.  An  attempt  was  made  at  this 
through  the  use  of  a graphite-compacted  material,  held  together  by  means 
of  two  axles  and  a set  of  steel  boots  at  each  end  of  the  composite  struc- 
ture. Figures  4 and  5 illustrate  this  concept,  which  was  eventually 
fabricated  and  tested. 
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Figure  3.  ROVAC  Vane  Loading  and  Deflections 
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Although  the  above  concept  did  offer  some  improvement,  load/deflection 
problems  were  still  prevalent  and  were  severely  inhibiting  the  develop- 
ment of  an  efficient  air-cycle  machine.  Based  on  similar  problems  and 
experience  gained  in  the  area  of  advanced  composites  for  bearings  and  j 

landing-gear-component  design,  it  was  suggested  that  advanced  filamentary  i 

composite  concepts  be  considered  for  use  as  a vane  material.  It  is  this  i 

consideration  which  led  to  the  formulation  of  the  program  discussed  within  i 

this  report. 

This  program  is  directed  toward  the  design  and  analysis  of  various 
lighter  weight  vane  concepts  employing  advanced  composite  materials. 

The  analysis  methods  used  employed  finite  element  techniques  and  refined 
modeling  methods.  Basically,  the  program  involved  five  primary  objectives 
which  were  as  follows; 

A)  Develop  an  automated  finite  element  vane  modeling  (AVM)  program 
from  which  various  vane  types,  sizes,  and  configurations  can  be  modeled. 

B)  Run  several  baseline  checks  through  NASTRAN  or  models  generated 
by  means  of  the  AVM  program.  These  models  will  consist  of  the  base 
graphite-compact  vane  (Figures  4 and  5)  to  which  analysis  results  can  be 
compared  to  available  analytical  and  experimental  results  on  the  vane 
itself. 

C)  Run  eight  direct  composite-substitution  analysis  surveys  to 
determine  the  effects  of  using  advanced  composites  within  the  vane  in  its 
current  design  state. 

D)  Run  several  additional  composite  runs  beyond  the  current  design 
state  whereby  the  negative  effects  of  incorporating  steel  boots  and  axles 
might  be  eliminated. 

E)  Based  on  the  results  of  C and  D above,  conduct  a series  of  runs 
incorporating  vane  thickness  variations  on  the  two  most  promising  designs. 
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SECTION  II 

NASTRAN  MODELING  PROGRAM 

Based  on  previous  efforts  in  the  area,  the  requirement  for  some  form 
of  aut.  ’’ted  data  generator  became  quite  obvious.  Two  primary  reasons 
existf  ich  determined  this  requirement  and  the  approach  to  be  selected 
for  tFe  «Ui.omated  data-generation  system.  The  first  was  the  fact  that 
the  sizing  of  the  air-cycle  machine  itself  involves  developmental  changes 
along  with  changes  in  vane  dimensions  and  sizes.  As  a result,  by  the  time 
one  vane  could  be  set  up  manually,  the  data  itself  became  obsolete.  The 
second  and  most  important  reason  was  the  fact  that,  through  a manual  mode 
of  data  generation,  little  other  than  design  verification  studies  could  be 
made  within  the  time  frames  involved.  Studies  of  this  type  were  very 
beneficial  but  what  was  really  required  were  parametric  analysis  surveys 
where  design  changes  could  be  studied  and  actually  implemented. 

The  actual  parametric  data-generation  system,  herein  referred  to  as  the 
AVM  (Automated  Vane  Modeling)  program,  employed  in  this  effort  involves  a 
total  of  eight  parameters  which  can  be  varied  in  the  automated  mode.  These 
are  as  follows: 

1.  The  ability  to  responsively  change  the  basic  dimensions  of  the 
vane  itself. 

2.  The  ability  to  vary  the  level  of  detail  of  the  analysis  effort 
(course  versus  refined  finite-element  models). 

3.  The  ability  to  incorporate  materials  changes  for  the  basic  vane 
itself. 

4.  The  ability  to  introduce  materials  variations  (multi-materials) 
within  the  design  of  a single  vane. 

5.  The  ability  to  include  the  use  of  advanced  composites  for 
reductions  in  vane  weight. 

6.  The  ability  to  vary  the  method  of  bearing  attachment  to  the  vane. 
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7.  The  ability  to  vary  the  method  of  load  transmission  through  the 
vane  by  either  double  axles,  a single  axle,  or  only  a partial 
axle. 

8.  The  ability  to  automatically  include  and  vary  all  of  the  symmetric 
boundary  constraints,  rotational  force  vectors,  panel  inertia 
variables,  and  vane  loading  and  reaction  boundary  conditions. 

1.  BASIC  MODEL 

The  AVM  program  developed  for  this  effort  consists  of  generating  the 
required  NASTRAN  input  for  seven  basic  data  sets.  Specifically,  these 
include  the  "GRID"  set,  "CQUAD2"  set,  "SPCl"  and  "SPCADD"  set  involving 
both  "RFORCE"  (rotational  forces),  "GRAV"  (Inertia  forces)  and  miscel- 
laneous external  forces  to  the  vane. 

The  basic  model  itself  consists  of  a rectangular  area  (Figure  6) 
divided  into  a set  of  NPX  (x  projection)  and  NPY  (y  projection)  lines. 
Typically,  this  results  in  a model  involving  a total  of  (NPX)  (NPY)  nodes 
in  the  structural  analysis  model. 

To  meet  the  eight  requirements  of  the  study,  the  uniform  distribution 
noted  in  Figure  6 was  further  developed  to  a point  where  the  actual  x and 
y values  for  each  x or  y projection  are  individually  input.  A typical 
resulting  model  based  on  this  setup  is  further  illustrated  in  Figure  7 
along  with  the  element  idealization. 

2.  "GRID"  DATA  SET 

In  accordance  with  the  required  NASTRAN  input  format,  the  following 
data  is  included  for  each  grid  point  within  the  model. 


ITEM  NO. 

1 

2 

3 

1 

4 

5 

6 

7 

8 

cor.  LOC, 

1-8 

9-16 

25-32 

33-40 

41-48 

1 

49-56 

57-64 

VARIABLE 

mm 

CP 

X 

y 

PS 

10 

L - ^ ^ 


NPY 
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I 

' Item 

Item 

Item 

Item 

Item 

j Item 

j Item 

! I tern 


1 - Alpha  Name  Grid 
Z - Grid  Point  Number  (integer) 

3 - Blank  (for  ROVAC  case) 

4 1 

cl  X,  y,  z values  for  Node  Point 
g)  (Floating  Point) 

g|  Blank  Fields  (for  ROVAC  case) 


3.  "CQUAD2"  DATA  SET 


It  is  within  the  "CQUAD2"  set  where  the  material  properties,  property 
variations  (multi-materials  per  vane),  and  the  use  of  advanced  composites 
are  introduced  into  the  vane  model.  Specifically,  these  parameters  are 
noted  as  variables  PID  and  TH  on  the  following  NASTRAN  data  set. 


ITEM  NO. 

1 

2 

3 

1 

4 

5 

6 

7 

8 

COL.  LOC. 

1-8 

m 

17-24 

25-32 

33-40 

41-48 

49-56 

57-64 

VARIABLE 

CQUAD2I 

1 EID 

PID 

G1 

G2 

G3 

G4 

TH 

Alphanumeric  name  "CQUAD2" 

Element  identification  no.  (EID)  (integer) 

Property  identification  no.  (PID)  (integer) 

(See  following  PID  discussion) 

Nodal  or  grid  point  connection  data  for  element  no.  EID 
Numbering  starts  at  lower  left  and  proceeds  counterclockwise 
around  element  (integers) 

Material  property  orientation  angle  for  advanced  composite 
type  materials  (see  following  TH  discussion) 

4.  PROPLRrV  IDENTIFICATION  DETERMINATION  (PID) 

Within  a typical  ROVAC  vane,  various  properties  exist  across  the 
vane  surface  in  the  form  of  end  plates,  boots,  adaptors,  or  varying 
materials  within  the  structure  itself.  To  accommodate  these  changes. 


Item  1 - 
Item  2 - 
Item  3 - 

i I tern  4 

j I tern  5 

! I tern  6 

! I tern  7 

I Item  8 - 
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the  AVM  program  was  set  up  to  include  four  primary  property  areas.  Typi- 
cally, these  property  areas  are  noted  in  Figure  8 as  PID  values  1 through  4, 

where: 

X.p,  Xnp,  X_p,  Xr,p,  and  Xcp=  input  values  (floating  point) 

^ ^ (FORTRAN  input  = XPl  through  XP5) 

All  CQUADS  lying  between  X^p  and  Xgp  will  have  PID  = 1 

All  CQUADS  lying  between  Xgp  and  X^p  will  have  PID  = 2 

All  CQUADS  lying  between  X^p  and  X^p  will  have  PID  = 3 

All  CQUADS  lying  between  X^p  and  X^p  will  have  PID  = 4 

5.  THETA  DETERMINATION  (TH) 

One  very  promising  method  of  significantly  reducing  the  vane  weight 
within  an  air-cycle  system  is  through  the  use  of  advanced  composite  mate- 
rials. To  include  the  possibility  for  considering  materials  of  this  type, 
a total  of  four  theta  regions  is  included  for  automated  setup  on  the  ROVAC 
vane  model.  Typically,  these  regions  are  shown  in  Figure  9 as  TH  (or 
theta)  values  THl  ...  TH4, 

where: 

X^j^ , ^CTH * ^DTH*  ^ETH  ” input  values  (floating  point) 

THl,  TH2,  TH3,  TH4  = input  values  (floating  point) 

(FORTRAN  input  = XTHl  through  XTH5) 

( and  THl  through  TH4) 

All  CQUADS  lying  between  and  Xgy^  will  have  TH  = THl 

Xgj^  and  Xgjg  will  have  TH  = TH2 

Xqj^  and  Xpj^  will  have  TH  = TH3 

Xg.j.g  and  X^y^  will  have  TH  = TH4 

where  THl  ...  TH4  = 0,  0 . 

1 4 

= material  property  orientation  angles  defining  at  which 
angle  the  anisotropic  or  orthotropic  properties  are  related. 
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6.  "CBAR"  DATA  SET  (EDGE  OPTION) 

The  mathematical  incorporation  of  double,  single,  or  partial  axle 
systems  through  the  vane  center  has  been  effected  by  the  use  of  NASTRAN 
bar  type  elements.  Typically,  a single  or  double  row  of  "CBAR"  type 
elements  extends  across  the  face  of  the  vane  model,  one  of  which  attaches 
to  a vane  edge  closeout  connected  to  a bearing  inner  race  structure. 

Figure  10  illustrates  segregated  bar  elements  for  a two-axle  structure 
and  Figure  11  shows  the  same  bar  elements  as  they  exist  in  the  total  model. 
Specific  locations  of  the  height  of  each  axle  is  defined  by  the  parameters 
YRODl  and  YR0D2.  A third  parameter  XRODl  is  used  to  define  the  location 
of  the  bearing  inner  race  edge.  Output  of  the  AVM  program  consists  of 
four  bar  elements  comprising  the  bearing  inner  race,  three  bar  elements 
comprising  the  attachment  to  the  vane  closeout,  six  bar  elements  simulating 
the  vane  closeout  stiffness,  and  a set  of  n-bar  elements  representing  the 
upper  and  lower  axles.  In  the  event  that  either  of  the  axle  y positions  do 
not  fall  coincident  to  an  existing  line  of  nodes  (i.e.,  nodes  6,  7,  8,  9, 

10,  Figure  11)  the  AVM  program  will  automatical ‘'y  generate  a set  of  nodal 
points  to  which  bar  elements  may  be  connected.  It  should  be  noted,  however, 
that  very  high  aspect  ratio  elements  could  result  when  an  axle  position  is 
very  close  to  but  not  coincident  with  an  existing  line  of  nodes.  In 
accordance  with  the  required  NASTRAN  input  format  the  following  data  is 
included  for  each  "CBAR"  element  within  the  vane  model. 


ITEM  NO. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

COL.  LOC. 

1-8 

9-16 

17-24 

25-32 

33-40 

41-48 

49-56 

57-64 

65-72 

VARIABLE 

CBAR 

EID 

PID 

GA 

GB 

XI 

X2 

1 

X3 

F 

1 

Item  1 - Alpha  name  "CBAR" 

Item  2 - Element  number  (EID);  sequenced  starting  after  all 
"CQUAD"  elements  have  been  numbered 

Item  3 - Property  identification  number 
(see  following  PID  discussion) 
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Figure  10.  Segregated  AVM  Bar  Elements  For  Double  Axle 
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37 

31 

38 

32 

39 

40 

35 

30 

25 

31 

26 

32 

27 

33 

28 

34 

^49 

26 

22 

.27 

23 

28 

24 

29 

47  48  49 

“A 

38 

,^48 

21 

39 

18 

22 

40 

19 

23 

4'l 

20 

24 

i 

25 

20 

YR« 

15 

YR0D2 
.0  1 

k 

)DI 

L_ 

47— 

13 

.16 

14 

17 

15 

18 

16 

19 

46  — 

9 

1 1 

10 

12 

1 1 

13 

12 

14 

^ vaniM 

5 

6 

7 

7 

8 

8 

9 ^^5 

1 ' 

2 ^ 

3 ^ 

4 ^ 

5 ' 

1 — ^ “j; 

Figure  11.  AVM  Bar  Elements  Within  Total  Model 
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Item  4 1 Nodal  or  grid  point  connection  data  for  element.  Numbering 
Item  Sf  proceeds  left  to  right  for  all  horizontal  elements  and 
* from  bottom  to  top  on  all  non-horizontal  elements. 


Item  6 
Item  7 
Item  8 


Vector  component  definition  for  bar  element  coordinate 
system  orientation. 

} = 


(0. 

for 

horizontal  bars) 

(1. 

for 

non-horizontal  bars) 

(1. 

for 

horizontal  bars) 

(0. 

for 

non-horizontal  bars) 

(0. 

for 

all  bars) 

} = 


= X2 
= X3 


Item  9 F=1  = flag  denoting  that  XI,  X2,  and  X3  will  be  read  as 

above 


7.  DOUBLE,  SINGLE,  OR  PARTIAL  AXLE  OPTION  (PID  DETERMINATION) 

Within  the  analysis  model,  provision  has  been  made  to  generate  either 
a single  or  double  axle  model.  For  the  single  axle,  the  upper  bearing 
attachment  axle  will  always  be  generated  without  the  lower  axle.  Each  of 
the  two  axles  in  the  two-axle  case  have  been  set  up  independently  of  each 
other  with  four  sets  of  material  properties  as  noted  in  the  PID  field 
(Item  3).  Soecifically,  within  the  model  the  upper  axle  is  defined  with 
PID's  of  5,  6,  and  7 where  the  lower  axle  is  defined  with  a PID  equal  to  8. 
Within  the  immediate  vane  region  constant  properties  have  been  assumed 
across  the  entire  axle.  Out  of  the  direct  vane  region,  the  closeout  edge 
has  been  assigned  a PID  equal  to  6 while  the  bearing  support  region  consists 
of  a PID  equal  to  5.  One  "CBAR"  element  in  the  bearing  structure  but 
immediately  adjacent  to  the  vane  upper  axle  also  carries  a PID  value  equal 
to  7,  in  that  it  is  upper  axle  material.  A typical  "CBAR"  set  for  the 
model  illustrated  in  Figures  10  and  11  is  noted  in  Table  1. 

For  a partial  or  no  axle  model  either  entire  "CBAR"  card  sets  can  be 
extracted  from  the  resulting  output  or  property  values  revised  to  negate 
the  effects  of  the  "CBAR"  elements. 
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TABLE  1 

"CBAR"  OUTPUT  FOR  FIGURES  10  AND  11 


CBAR 

33 

5 

46 

47 

0.00 

1.00 

0.00 

1 

CBAR 

34 

5 

47 

48 

0.00 

1.00 

0.00 

1 

CBAR 

35 

5 

48 

49 

0.00 

1.00 

0.00 

1 

CBAR 

36 

5 

49 

50 

0.00 

1.00 

0.00 

1 

CBAR 

37 

7 

50 

26 

0.00 

1.00 

0.00 

1 

CBAR 

38 

7 

26 

27 

0.00 

1.00 

0.00 

1 

CBAR 

39 

7 

27 

28 

0.00 

1.00 

0.00 

1 

CBAR 

40 

7 

28 

29 

0.00 

1.00 

0.00 

1 

CBAR 

41 

7 

29 

30 

0.00 

1.00 

0.00 

1 

CBAR 

42 

8 

6 

7 

0.00 

1.00 

0.00 

1 

CBAR 

43 

8 

7 

8 

0.00 

1.00 

0.00 

1 

CBAR 

44 

8 

8 

9 

0.00 

1.00 

0.00 

1 

CBAR 

45 

8 

9 

10 

0.00 

1.00 

0.00 

1 

CBAR 

46 

6 

11 

16 

1.00 

0.00 

0.00 

1 

CBAR 

47 

6 

16 

21 

1.00 

0.00 

0.00 

1 

CBAR 

48 

6 

21 

26 

1.00 

0.00 

0.00 

1 

CBAR 

49 

6 

26 

31 

1.00 

0.00 

0.00 

1 

CBAR 

50 

6 

31 

36 

1.00 

0.00 

0.00 

1 

CBAR 

51 

6 

36 

41 

1.00 

0.00 

0.00 

1 

CBAR 

52 

5 

21 

50 

1.00 

0.00 

0.00 

1 

CBAR 

53 

5 

50 

31 

1.00 

0.00 

0.00 

1 

8.  "CBAR"  SET  (OPTION  NON-EDGE) 

In  the  previous  axle  model,  an  existing  set  of  nodes  along  an  existing 
y edge  was  used  to  define  the  axle  bar  elements  or  one  was  automatically 
generated.  For  some  cases,  however,  it  may  be  required  that  the  axle  be 
modeled  across  the  centroid  or  some  portion  of  the  element  other  than  the 
element  edge.  For  the  latter  case,  an  entirely  new  set  of  nodal  values 
is  generated.  Within  the  AVM  program  nodes  may  be  suppressed  and  the 
generation  of  a line  of  non-edge  nodes  can  be  interchanged.  Figure  12a 
illustrates  this  suppressed  mode  for  a single  axle  non-edge  generation 
scheme.  Node  points  51  through  55  have  been  added  to  the  base  model  and 
from  which  the  axle  noted  in  Figure  12b  was  modeled.  All  of  these  nodes 
are  totally  independent  of  any  element  edge  and  truly  lie  across  one  of 
the  existing  elements  in  the  structural  model.  It  becomes  obvious, 
however,  that  the  axles  modeled  in  this  case  are  not  inherent  to  the  basic 
structural  model.  Some  additional  modeling  will  therefore  be  required  to 
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relate  the  axles  to  the  model.  This  additional  modeling  has  been  left  up 
to  the  user  and  has  not  been  automated  into  the  AVM  program. 

For  the  non-edge  case,  the  vane  closeout  edge  is  modeled  such  that 
two  bars  are  modeled  across  the  left  edge  of  the  element  divided  by  the 
axle.  Two  additional  bars  are  generated  above  and  below  the  inner  race 
attachment  points  as  is  done  in  the  edge  case.  All  other  modeling  para- 
meters are  identical  to  the  edge  case. 

9.  "SPC"  DATA  SET  (MATHEMATICAL  CONSTRAINTS) 

Within  the  AVM  program  various  SPC  (Single  Point  Constraint)  sets  are 
also  generated.  These  sets  include  reaction  points,  unidirectional  bearing 
constraints,  constraints,  and  symmetry  constraints  for  symmetrical 
modeling.  Each  of  these  SPC  sets  is  discussed  in  the  following  four 
sections. 

a.  Symmetry  Condition 

Due  to  the  fact  that  the  vane  is  symmetrically  loaded  and 
dimensionally  identical  about  its  centerline,  only  one  half  of  the  vane 
need  be  modeled.  To  effect  this  symmetry  the  AVM  program  selects  the 
appropriate  edge  nodes  along  the  centerline  (Figure  13)  and  generates  an 
appropriate  NASTRAN  "SPCl"  card  for  each  node  point  along  this  centerline 
edge.  The  actual  constraints  used  resist  both  motion  in  the  x direction 
and  rotation  about  the  y axis.  The  specific  SPCl  format  output  is  as 
follows: 


ITEM  No. 

1 

2 

3 

4 

5 

COL.  LOC. 

1-8 

9-16 

17-24 

25-32 

33  40 

VARIABLE 

SPCl 

SID 

C 

G 

I 
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Item  1 - Alphanumeric  name  "SPCl" 

Item  2 - Constraint  ID  number,  initialized  at  SID  = 100  (integer) 
followed  by  consecutive  +1  incrementing 

Item  3 - Constraint  factor  no.  15  denoting  a constraint  in  the 
1 and  5 or  X and  directions,  respectively  (integer) 

Item  4 - Nodal  point  to  be  constrained  (integer) 

Item  5 - N/A  (blank) 


I — I 


Figure  13.  AVM  Symmetrical  SPC  Set  Definition 
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b.  Unidirectional  Bearing  Constraint 

Within  the  analysis  model  a "SPCl"  constraint  is  generated  repre- 
senting the  total  restraint  existing  (in  the  y plane)  at  the  bearing.  At 
this  point  y motion  is  reacted  by  the  ROVAC  track,  x motion  is  reacted  by 

the  machine  end  plates,  M , and  M are  reacted  through  the  rotor  and  z is 

X y 

fixed  because  of  the  y plane  constraint.  Typically,  this  bearing  con- 
straint is  always  applied  to  the  left  most  nodal  point  (No.  37,  Figure  13) 
on  the  inner  race  axle  model.  For  this  case,  one  "SPCl"  card  is  generated 
with  C = 123456.  The  assigned  SPC  set  number  for  this  card  is  500. 

c.  Reaction  Points 

Primary  loads  within  the  ROVAC  vane  structure  are  generated  from 
both  rotational  and  linear  inertia  effects.  Typically,  these  loads  are 
reacted  in  the  ^ direction  at  the  bearing  inner  race  mount.  This  mount 
in  the  automated  generator  consists  of  the  four  left  most  nodal  points 
(i.e.,  37  through  40,  Figure  13).  To  satisfy  the  ^ reaction,  three  of 
these  four  points  must  be  additionally  constrained.  This  is  accomplished 
through  the  use  of  a "SPCl"  card  containing  all  three  points  where  C is 
set  equal  to  2.  The  SPC  set  number  for  this  card  is  502. 

d.  Z Axis  Moment  Constraint  (M^) 

The  standard  "CQUAD2"  NASTRAN  constraint  set  has  been 
represented  in  the  vane  model  by  the  SPC  set  number  501.  Within  this  set 
all  nodal  points  within  the  model  are  constrained  in  the  R3  direction. 

A point  of  caution  or  possible  error  is  presented  at  this  point,  in  that 
those  "CBAR"  elements  requiring  M^  freedom  could  adversely  be  affected 
through  this  SPCl  imposition.  If  this  is  the  case,  a separate  hand- 
generated SPCl  set  would  have  to  be  developed  for  which  would  not 
include  nodal  points  associated  with  any  "CBAR"  type  elements. 

e.  Combined  Constraints 

All  of  the  SPC  sets  generated  for  the  vane  model  are  combined 
into  a single  constraint  set.  This  constraint  set  is  assembled  by  means 
of  the  "SPCADD"  card  and  continuation  cards.  The  final  resulting  SPC  set 
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for  NASTRAN  case  control  is  SID  = 99.  Table  2 is  a tabulation  of  all  the 
SPC  sets  and  SPC  combination  for  the  model  noted  in  Figure  13. 

TABLE  2 

AVM  GENERATED  CONSTRAINT  SET 


SPCl 

100 

15 

4 

SPCl 

101 

15 

8 

SPCl 

102 

15 

12 

SPCl 

103 

15 

16 

SPCl 

104 

15 

20 

SPCl 

105 

15 

24 

SPCl 

106 

15 

28 

SPCl 

107 

15 

32 

SPCl 

108 

15 

36 

SPCl 

109 

15 

45 

SPCl 

500 

123456 

37 

SPCl 

501 

6 

1 

THRU 

45 

SPCl 

502 

2 

38 

39 

40 

SPCADD 

99 

100 

101 

102 

103 

104 

105 

106 

ABCIO 

+BC10 

107 

108 

109 

500 

501 

502 

(Reference  Figure  13) 

10.  PROPERTY  AND  MATERIALS  DATA 

Quadrilateral  properties  within  the  ROVAC  model  are  defined  by  means 
of  "PQLIAD2"  cards  representing  each  of  the  four  vane  areas  previously 
discussed. 

a.  "PQUAD2"  Data 

The  specific  "PQUAD2"  data  output  for  the  ROVAC  vane  model  is  as 

follows: 


ITEM  NO. 

1 

2 

3 

4 

COL.  LOC. 

1-8 

9-16 

17-24 

25-32 

VARIABLE 

PQUAD2 

PID 

MID 

T 

Item  1 - Alphanumeric  name  "PQUAD2" 

Item  2 - Property  identification  no.  (PID)  (INTEGER) 

Item  3 - Material  identification  no.  (references  a specific  material 
card  where  the  specific  material  properties  are  defined  for 
the  property  identification  (PID)  (INTEGER) 

Item  4 - Thickness  of  all  elements  relating  to  this  property  (Real  Number) 
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b.  "PBAR"  Data 

Bar  element  properties  within  the  ROVAC  model  are  defined  by 
means  of  four  "PBAR"  cards  representing  each  of  the  four  axle  or  bearing- 
support  regions  previously  discussed.  Due  to  the  fact  that  some  of  the 
properties  are  small  and  are  self-computed  by  the  automated  data-generator 
program,  a la’^ge  field  card  was  required.  The  specific  large  field  "PBAR" 
data  output  for  the  ROVAC  vane  model  is  on  two  cards  as  follows: 


Item  1 - Alphanumeric  name  "PBAR*"  where  the  * flags  a large  field 
card  and  the  following  four  fields  of  data  on  the  card  will 
be  read  in  fields  of  16  rather  than  8. 

Item  2 - Bar  property  identification  number  (PID)  (Integer) 

Item  3 - Material  ID  No.  (Ref.  PQUAD2  card) 

Item  4 - Bar  cross-sectional  area  (Real) 

Item  5 - Area  moment  of  inertia  (Direction  1)  (Real) 

Item  6 - Coded  continuation  card  flag 
Item  7 - Continuation  cards  related  to  coded  item  number  6 
(*  = large  field  flag  for  items  8 and  9) 

Item  8 - Area  moment  of  inertia  (Direction  2)  (Real) 

Item  9 - Torsional  constant  (Real) 
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c.  "MAT!"  Data 

Each  of  the  "PQUAD2"  and  "PBAR"  cards  noted  previously  contains 
a material  identification  number  as  Item  3.  For  the  ROVAC  vane  each  of 
these  ID  numbers  is  unique  for  each  of  the  eight  property  cards.  As  a 
result,  eight  unique  materials  cards  are  required  for  the  ROVAC  model. 
Typically,  each  of  these  eight  materials  references  are  output  as  follows: 


ITEM  NO. 

1 

3 

4 

5 

6 

7 

COL.  LOC. 

1-8 

9-16 

17-24 

25-32 

33-40 

41-48 

49-56 

VARIABLE 

E 

G 

■■ 

Item  1 - Alphanumeric  name  "MATl" 

Item  2 - Materials  property  identification  number  per  property 
card  reference  (Integer) 

Item  3 - Young's  modulus  (Real) 

Item  4 - Blank  for  ROVAC  Case  (computed  within  NASTRAN) 

Item  5 - Poisson's  ratio  (Real) 

Item  6 - Mass  density  (Real) 

Item  7 - Blank  for  ROVAC  model 

11.  ROVAC  VANE  LOADING 

As  stated  previously,  vane  loads  are  generated  from  two  primary 
sources,  consisting  of  rotational  and  inertia  forces  generated  during 
air-cycle  machine  operation.  Specifically,  the  rotational  forces  result 
from  pure  rotor  rotation,  whereas  the  inertial  force  is  added  due  to  the 
effects  of  in  and  out  movement  of  the  vane  through  the  rotation  angle. 

To  simulate  these  loading  values  in  the  ROVAC  model,  two  distinct  input 
sets  are  generated  by  the  AVM  program. 

a.  Rotational  Force  ("RFORCE") 

The  first  load  condition  consists  of  a rotation  force  ("RFORCE") 
data  set  defining  a static  loading  condition  due  to  a centrifugal  force 
field.  Based  on  the  mass  properties  of  the  model,  a corresponding 
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centrifugal  loading  for  each  element  in  the  model  will  be  generated.  The 
specific  ROVAC  AVM  output  for  the  "RFORCE"  data  is  as  follows: 


ITEM  NO. 

1 

2 

3 

4 

5 

6 

COL.  LOC. 

1-8 

9-24 

25-40 

41-56 

57-72 

73-80 

VARIABLE 

RFORCE* 

SID 

G 

CIO 

A 

QRFORCEl 

ITEM  NO. 

7 

8 

9 

10 

COL.  LOC. 

1-8 

9-16 

17-24 

25-32 

VARIABLE 

+RF0RCE1 

N1 

Item  1 - 


Item  2 - 

Item  3 - 

Item  4 - 

Item  5 - 

Item  6 - 
Item  7 - 
Item  8 
Item  9 
Item  10 


Alpha  name  "RFORCE*"  where  the  * flags  a large  field  set 
and  the  following  4 fields  of  data  will  be  read  in  fields 
of  16  rather  than  8. 

Unique  load  set  identification  number  (Integer)  (for 
ROVAC  SID  = 20). 

Grid  point  identification  number  indicating  a center  of 
rotation  for  the  ROVAC  vane  model  (Integer). 

Blank  for  ROVAC  case  indicating  that  the  basic  co-ordinate 
system  defines  the  rotation  direction. 

Scale  factor  for  rotational  velocity  in  revolutions  per 
unit  time  (ROVAC  = rps). 

Coded  continuation  card  flag. 

Continuation  card  related  to  coded  Item  No.  6. 

Rectangular  components  of  rotation  direction  vector  (Real). 


b.  Inertia  Loading  ("GRAV") 

The  second  load  condition  consists  of  a gravity  vector  "GRAV" 
data  set  defining  a gravity  vector  simulating  an  inertia  type  loading  for 
the  ROVAC  vane  model.  Based  on  the  mass  properties,  then,  a corresponding 
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inertia  loading  for  each  element  in  the  model  will  be  generated.  The 


specific  ROVAC 

AVM  output  for  this  "GRAV"  data  is  as  follows: 

ITEM  NO. 

1 

2 

3 

4 

5 

6 

COL.  LOC. 

1-8 

9-24 

25-40 

41-56 

57-72 

73-80 

VARIABLE 

GRAV* 

SID 

CID 

G 

N1 

QGRAVl 

ITEM 

NO. 

7 

8 

9 

COL. 

LOC. 

1-8 

9-16 

17-24 

VARIABLE 

+GRAV1 

N2 

N3 

Item  1 - Alpha  name  "GRAV*"  where  the  * flags  the  large  field  set 
noted  previously. 

Item  2 - Unique  load  set  identification  number  (Integer) 

(for  ROVAC  Case  SID  = 21). 

Item  3 - Coordinate  system  identification  number  set  to  blank  (0) 

for  the  ROVAC  case.  (References  the  basic  coordinate  system.) 

Item  4 - Actual  gravity  vector  value  (Reference  the  following 
gravity  vector  discussion). 


Item  5 
Item  8 - 
Item  9 
Item  6 - 
Item  7 - 


Gravity  vector  components  (Real). 

Coded  continuation  card  flag. 

Continuation  card  related  to  coded  item  No.  6. 


c.  Gravity  Vector  (Item  4) 


Within  the  air-cycle  machine,  it  was  noted 
traveling  about  an  elliptical  path,  which  results  in 
lation  of  the  vane  itself  throughout  each  rotational 
and  direction  of  this  travel  is  of  course  controlled 
track  and  the  slots  within  the  rotor  are  as  noted  in 
1 and  2.  As  a result  of  this  in-and-out  vane  motion 


that  the  vane  is 
an  in-and-out  trans- 
cycle. The  magnitude 
by  the  elliptical  cam 
Section  I,  Figures 
at  high  rotor  speeds. 


I 


28 


AFFDL-TR-77-49 


a very  significant  inertia  load  is  generated  on  the  vane.  In  addition 
to  this  load,  a friction  load  is  further  imposed  as  a result  of  the  vane 
sliding  in-and-out  of  the  rotor  slot.  Both  of  these  forces  are  repre- 
sented within  the  AVM  NASTRAN  model  by  means  of  an  acceleration  or  gravity 
vector  value  on  a GRAV  input  card.  The  actual  value  of  this  vector  can 
be  derived  from  several  sources  or  it  could  be  computed  from  the  machine 
geometrical  parameters.  The  provision  for  computing  the  value  is  provided 
by  means  of  eight  input  parameters  ALDATA(l)  through  ALDATA(8).  For  this 
provision,  the  specific  geometric  equation  would  have  to  be  added  to  the 
program.  For  the  analysis  conducted  in  this  report,  a representative 
input  value  was  derived  from  an  additional  computer  program  which  specifi- 
cally analyzes  the  geometric  relationship.  The  actual  AVM  input  value  is 
an  acceleration  value  noted  in  inches  per  second. 

d.  Combined  Loads 

Both  the  "RFORCE"  load  and  "GRAV"  load  sets  are  combined  into  a 
common  load  set  through  the  use  of  a simple  load  card.  For  the  ROVAC 
case,  the  rotational  and  gravity  force  sets  have  been  identified  as  20  and 
21,  respectively.  The  resulting  combined  load  set  has  been  identified  as 
number  30  for  NASTRAN  case  control.  The  specific  ROVAC  AVM  output  for  the 
"LOAD"  set  is  as  follows; 


ITEM  NO.. 

1 

2 

3 

4 

5 

6 

7 

COL.  LOC. 

1-8 

9-16 

17-24 

25-32 

33-40 

41-48 

49-56 

VARIABLE 

LOAD 

SID 

S 

SI 

LI 

L2 

Item  1 - Alpha  name  "LOAD" 

Item  2 - Load  set  ID  number  (SID  - 30  for  ROVAC  AVM) 

Item  3 - Scale  factor  for  Item  No.  2 load  set  (S  = 1.000  for  ROVAC  AVM) 

Item  5»-  Load  set  identification  numbers  to  be  combined  into  the  SID 
Item  7l-  set.  For  the  ROVAC  case,  LI  = 20  and  L2  = 21 . 

Item  4i-  Respective  scale  factors  (SI,  S2)  for  load  sets  (LI,  L2). 

Item  6^-  for  ROVAC  AVM  SI  = S2  = 1.000. 
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e.  Parameter  Usage  ("WTMASS"  and  "GRDPNT") 

In  accordance  with  NASTRAN  requirements  and  desired  ROVAC  output 
two  "PARAM”  cards  are  generated  by  the  AVM  program.  These  include  a grid 
point  "GRDPNT"  and  a weight  mass  "WTMASS"  parameter  in  accordance  with 
the  following; 


ITEM  NO. 

1 

2 

3 

COL.  LOC. 

1-8 

9-16 

17-24 

CARD  1 

VARIABLE 

PARAM 

WTMASS 

.002588 

CARD  2 

VARIABLE 

PARAM 

GRDPNT 

0 

The  weight  mass  parameter  ("WTMASS")  will  result  in  all  structural 
mass  matrix  terms  being  multiplied  by  the  value  contained  in  Item  3.  For 
the  ROVAC  case,  the  mass  parameter  is  1/g,  resulting  in; 


where 

1/g  = 1/(32. 2)(12)  = .002588  weight  mass  parameter 
F = force 

(jj  = angular  velocity 
r = radius 

The  grid  point  ("GRDPNT")  parameter  will  result  in  the  execution  of 
the  weight  generator  package  where  principal  vane  masses,  centers  of 
gravity,  and  inertias  about  the  center  of  gravity  will  be  generated. 

The  value  contained  as  Item  3 indicates  the  grid  point  to  be  used  as  a 
reference  point  for  this  information.  For  the  AVM  program,  this  value  has 
been  set  at  zero,  and  as  a result  the  reference  point  is  taken  as  the 
origin  of  the  basic  coordinate  system. 

f.  Rotational  Axis  Definition 

The  generation  of  rotational  forces  for  ROVAC  studies  require 
the  definition  of  a rotational  axis  about  which  the  vane  shall  rotate. 
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As  noted  previously  in  the  RFORCE  output.  Item  3,  (G)  is  required  to 
define  this  axis  specifically.  The  value  G is  a designated  grid  point 
lying  on  the  axis  of  rotation  from  which  the  centrifugal  force  is  derived. 
For  the  ROVAC  AVM  case,  this  grid  point  is  an  entirely  independent  point 
where: 


X = z = 0.0 
y = y^  = AXIS 

AXIS  = distance  from  vane  edge  to  the  center  of  rotation 
g.  Program  Termination 

For  program  termination,  the  required  "ENDOATA"  card  is  self- 
generated within  the  ROVAC  AVM  program. 
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SECTION  III 

AVM  PROGRAM  INPUT  AND  OPERATION 

As  noted  in  the  program  listing  contained  in  Appendix  A,  the  AVM 
program  consists  of  a main  program  called  VANE  and  seven  associated 
subroutines.  These  subroutines  and  their  respective  functions  are  as 
follows : 

READIN  - Program  to  simply  read  in  and  store  all  of  the  required 
program  input  data. 

GRID  - Program  to  generate  all  of  the  X-Y  co-ordinate  values 
associated  with  each  nodal  point  within  the  structural 
model,  (output  to  TAPE  7) 

CQUAD2  - Program  to  calculate  and  set  up  all  of  the  quadralateral 
elements  within  the  finite  element  model,  (output  to 
TAPE  7) 

LISTIN  - Proqram  to  write  out  all  input  data  and  associated  internally 
computed  input  values,  (output  to  TAPE  6) 

CBAR  - Program  to  calculate  and  set  up  all  of  the  bar  elements 
within  the  finite  element  model,  (output  to  TAPE  7) 

CALCIN  - Program  to  calculate  and  set  up  all  of  the  PQUAD2,  PBAR, 

MATI,  SPCl  (500),  SPCl  (501),  SPCl  (502),  RFORCE,  GRAV, 

PARAM  (WTMASS),  PARAM  (GRDPNT),  GRID  (Rotational),  and 
ENDDATA  cards  required  for  the  finite  element  analysis, 
(output  to  TAPE  7) 

SPC  - Program  to  calculate  and  set  up  all  of  the  constraint  cards 
required  to  model  the  symmetry  conditions  associated  with 
the  vane  model,  (output  to  TAPE  7) 

a.  Input  Subscripts 

AVM  program  input  variables  are  generally  related  to  eight 
sections  or  materials  property  areas  contained  within  the  model.  Four  of 
these  areas  are  contained  within  the  vane  itself  while  the  remaining  four 
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are  restricted  to  the  vane  axles.  These  eight  areas  are  related  to 
property  variables  within  the  program  by  means  of  subscripts  or  arrays 
numerically  equivalent  to  the  area  number.  Each  of  these  eight  areas 
is  noted  in  Figure  14  and  defined  as  follows: 

A1  through  A4  - Four  areas  across  vane  face  (left  to  right  where 
thickness,  p,  p,  E,  and  other  properties  may  be 
varied) . 

A5  - Upper  axle  bearing-support  area  extending  to  the 

left  of  the  axle  transition  area. 

A6  - Axle  transition  and  vane  closeout  areas  located  to 

the  immediate  left  of  the  vane  itself. 

A7  - Upper  axle  extending  through  the  vane  and  into  the 

upper  axle  bearing-support  area. 

A8  - Lower  axle  area  extending  through  the  vane, 

b.  Computed  Variables 

Within  the  AVM  program  axle  areas  and  axle  area  moments  of 

inertia  I and  I ) are  self  computed  and  automatically  programmed  into 
XX  yy 

the  NASTRAN  output  deck.  The  computations  are  based  entirely  on  AVM 
input  and  are  individually  related  to  axle  areas  5 through  8. 

Within  NASTRAN  vane  accelerations  and  inertia,  forces  are 
computed  based  on  a vane  speed  in  revolutions  per  second.  The  resulting 
NASTRAN  output  from  the  AVM  program  automatically  accounts  for  these 
units.  Gravity  acceleration  values  can  also  be  computed  within  the  AVM 
program  if  desired.  These  values  are  based  upon  air-cycle  machine  design 
criteria  as  defined  by  AVM  program  input  variables.  The  actual  compu- 
tation is  based  on  input  data  noted  as  ALDATA  (1)  through  ALDATA  (8). 
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c.  Program  Input 


Input  to  the  AVM  program  is  accomplished  by  means  of  card  or 
card  sets  containing  all  of  the  variables  required  for  AVM  program 
operation.  Each  of  the  individual  cards  or  card  sets  is  described  as 
follows : 


CARD  1 


CARD  2 
(SET) 


CARD  3 
(SET) 


CARD  4 


INPUT  VARIABLES  = NPX,  NPY,  PRINT 

FORMAT  = 215,  LI 

NPX  = number  of  points  to  be  input  along  the  X axis  of 
the  vane  from  which  specific  X values  will  be 
calculated  for  each  nodal  point. 

NPY  = same  as  NPX  only  for  Y axis  and  Y values 

PRINT  = T or  F value  indicating  whether  or  not  an  output 
listing  of  the  program  input  is  desired. 

INPUT  VARIABLES  = XIN(I)  where  I = 1 to  NPX 

FORMAT  = 8F10.0 

XIN  = actual  X values  to  be  input  for  the  X sectioning 
of  the  vane  structure  (values  in  inches). 

INPUT  VARIABLES  = YIN(I)  where  I = 1 to  NPY 

FORMAT  = 8F10.0 

YIN  = actual  Y values  to  be  input  for  the  Y sectioning 
of  the  vane  structure  (values  in  inches). 

INPUT  VARIABLES  = XPl , XP2,  XP3,  XP4,  XP5 

FORMAT  = 8 FI  0.0 

XPl  through  XP5  = 5 X values  which  divide  the  vane  into 
four  distinct  areas  for  the  future 
assignment  of  varying  material  properties 
across  the  vane  face  (values  in  inches). 
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CARD  5 - INPUT  VARIABLES  = XTHi , XTH2,  XTH3,  XTH4,  XTH5 

FORMAT  = 8F10.0 

XTHI'  through  XTH5  = 5 X values  which  divide  the  vane  into 

four  distinct  areas  for  the  future 
assignment  of  varying  theta  values 
defining  material  property  orientation 
angles  across  the  vane  face  (values 
in  inches). 

CARD  6 - INPUT  VARIABLE  = THl , TH2,  TH3,  TH4 

FORMAT  = 8F10.0 

THl  through  TH4  = actual  property  orientation  angle  (a) 

in  each  of  the  four  areas  defined  by 
card  number  4 (value  in  degrees  as 
measured  from  an  axis  along  the  bot- 
tom of  the  vane  from  a leftwise 
origin  where  up  is  a positive  value). 

CARD  7 - INPUT  VARIABLES  - NROD,  XRODl , YRODl , YR0D2,  NEOPTIN 

FORMAT  = 110,  3F10.0,  110 

NROD  = number  of  axles  to  be  modeled  within  the  vane 
structure  (restricted  to  1 or  2). 

XRODl  = left  most  origin  of  the  upper  axle  which  defines 
the  length  of  axle  extending  between  the  bearing 
and  vane  edge  (value  in  inches). 

YRODl  = distance  of  upper  axle  from  the  lower  edge  of 
the  vane  structure  (value  in  inches). 

YR0D2  = distance  of  the  lower  axle  from  the  lower  edge  of 
the  vane  structure  (value  in  inches). 

NEOPTIN  = defines  if  the  non-edge  option  is  to  be  selected 
(any  number  greater  than  or  equal  to  1 will  result 
in  a non-edge  model ) . 
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CARD  8 - 


CARD  9 - 

CARD  10  - 


CARD  11  - 


1 


INPUT  VARIABLES  - IFMT 

FORMAT  = 15 

IFMT  = flag  indicating  in  what  format  the  remaining  data 
is  to  be  read  in.  If  IFMT  is  not  equal  to  1 the 
program  will  skip  card  sets  9 through  13,  This 
option  is  for  the  input  of  more  detailed  properties 
in  the  format  required  for  composite  materials. 

The  option  currently  is  not  in  use  and  IFMT  must 
be  set  equal  to  one  (Value  = 1). 

INPUT  VARIABLE  = XJ5,  XJ6,  XJ7,  XJ8 

FORMAT  = 40X,  4F10.0 

XJ5  through  XJ8  = axle  torsional  constant  in  respective 
areas  5 through  8 as  noted  in  Figure  9 
(axle  properties  restricted  to  circular 
section) . 

INPUT  VARIABLES  = THICKl,  THICK2,  THICK3,  THICK4, 

RADIUS5,  RADIUS6,  RADIUS7,  RADIUS8 

FORMAT  = 8F10.0 

THICKl  through  THICK4  = vane  thickness  values  in  respective 

vane  areas  1 through  4 (values  in 
inches) 

RADIUS5  through  RADIUS8  = Radius  values  of  axle  components 

as  defined  by  the  axle  areas 
noted  in  Figure  9 (values  in  inches). 

INPUT  VARIABLES  - El,  E2,  E3,  E4,  E5,  E6,  E7,  E8 

FORMAT  = 8 (2X,  A8) 

El  through  E8  = modulus  (E)  values  for  vane  areas  1 
through  4 and  axle  areas  5 through  8. 

Input  is  in  E format  (i.e.,  3.0E7  for 
steel ) 
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CARD  12  - INPUT  VARIABLES  = XMUl , XMU2,  XMU3,  XMU4,  XMU5,  XMU6, 

XMU7,  XMU8 

FORMAT  = 8F10.0 

XMUl  through  XMU8  = Poisson's  Ratio  values  for  vane  areas 
1 through  4 and  axle  areas  5 through  8. 

CARD  13  - INPUT  VARIABLES  = RHOl , RH02,  RH03,  RH04,  RH05,  RH06, 

RH07,  RH08 

FORMAT  = 8F10.0 

RHOl  through  RH08  = density  values  for  vane  areas  1 through 
4 and  axle  areas  5 through  8 (values  in 
Ib/in^. 

CARD  14  - INPUT  VARIABLES  - VEL,  AXIS,  GRAV 
FORMAT  = E FI  0.0 

VEL  = compressor  vane  rotational  velocity  (value  in 
revolutions  per  minute). 

AXIS  = distance  from  vane  edge  to  center  of  rotation  for 
compressor  vane  (value  is  in  inches;  ref.  AXIS 
Figure  9). 

GRAV  = gravity  acceleration  value  input  at  the  option  of 

the  user.  If  GRAV  is  equal  to  zero,  a gravity 

acceleration  value  based  on  the  data  supplied  on 

2 

Card  No.  15  will  be  computed  (value  = in/sec  ). 

CARD  15  - INPUT  VARIABLES  = ALDATA(l),  ALDATA(2),  ALDATA(3),  ALDATA(4), 

ALDATA(5),  ALDATA(6),  ALDATA(7),  ALDATA(8)  (*) 

FORMAT  = 8F10.0 

ALDATA(l)  through  ALDATA(8)  = algorithm  data  variables  1 

through  8 for  computing  the 
gravity  acceleration  value  to 
simulate  the  vane  forces  due  to 
the  inertia  effects  resulting 
from  the  elliptical  path  or  rota- 
tion within  the  air  cycle  ma- 
chine. (*) 

♦Currently  not  set  up  within  the  AVM  program  (blank  card  required). 


38 


AFFDL-TR-77-49 


d.  Data  Summary 

A complete  summary  of  all  associated  input  variables  and  input 
locations  is  presented  in  Table  3. 

e.  Test  Cases 

Four  test  cases  were  set  up  to  check  out  the  capabilities  of 
the  AVM  program.  Each  of  these  cases  is  included  in  the  section  as 
sample  I/O  examples  of  the  AVM  program.  AVM  capabilities  in  addition  to 
those  noted  in  these  four  cases  should  be  individually  checked  out  by  the 
user  prior  to  running  within  NASTRAN. 

Two-axle  model  where  both  the  upper  and  lower 
axles  fall  on  existing  lines  of  nodes.  (Bearing- 
support  axle  at  bottom). 

Two-axle  model  where  the  upper  axle  falls  on 
an  existing  nodal  line  and  the  lower  axle 
requires  a self-generated  edge.  This  model  is 
also  representative  of  an  existing  vane  in  an 
air  cycle  machine. 

Two-axle  model  where  the  lower  axle  is 
designed  to  fall  between  two  nodal  lines. 

For  this  case,  a separate  axle  independent 
of  the  vane  structure  is  generated  resulting 
in  15  rows  of  elements  rather  than  the  16 
noted  in  Case  No.  2. 

Example  application  of  a non-edge  case  where 
a model  has  been  modified  to  include  3 cutout 
holes  in  the  vane  face. 

f.  Program  Input/Output 

Specific  input  data  for  each  of  the  test  cases  1 through  3 is 
noted  in  Table  4.  Test  case  4 is  not  included  because  this  case  simply 
extracted  "CQUAD"  and  "GRID"  cards  fi'orn  test  case  3 output  to  simulate 


CASE  1 Edge  Model 


CASE  2 Self- 
Generating  Edge 


CASE  3 Non-Edge 
Option 


CASE  4 Non-Edge 
Option  with  Cutout 
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the  lightning  holes  in  an  early  ROVAC  vane  design.  A final  version  of 
this  No.  4 model  with  triangular  elements  added  to  the  hole  area  is 
included  in  Appendix  B. 

Details  as  to  element  and  nodal  numbering  for  each  of  the  test  cases 
are  noted  in  the  twelve  calcomp  plots  included  in  Appendix  B.  Also  noted 
in  this  appendix  are  the  abbreviated  output  listings  for  each  of  the 
test  cases. 
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SECTION  IV 
VANE  ANALYSIS  SETUP 

Twenty-two  vanes  (total)  were  setup  within  the  program  and  analyzed  by 
means  of  both  NASTRAN  and  the  AVM  program.  The  overall  analysis  was 
conducted  in  three  parts,  which  ranged  from  a very  minimal  design  change 
to  a total  redesign  of  the  vane  itself.  The  first  part  consisted  of  a 
direct  substitution  analysis,  whereby  the  ROVAC  vane  structure,  with  the 
associated  bearing,  bearing  axle  supports,  and  vane  boot,  was  retained 
in  its  entirety.  The  second  and  third  parts  of  the  study  involved  the 
elimination  of  axles,  boots,  and  vane  thickness  restrictions  associated 
with  the  ROVAC  vane  design. 

1.  PART  I (ROVAC  AND  ADVANCED  COMPOSITES) 

The  Part  I analysis  consisted  of  ten  NASTRAN  models  representing  ten 
graphite  fiber  materials  variations  within  the  ROVAC  vane  design.  For 
this  design  the  concept  of  using  a boot  and  axle  system  (Figures  4 and  5) 
was  retained  and  varied  orthotropic  materials  properties  were  introduced 
which  effectively  replaced  the  existing  Graphite  Compact  vane.  The  ten 
materials  considered  within  the  analysis  survey  were  as  follows: 

CASE  1 - 


i 

i 

I CASE  2 - 


CASE  3 - 


1 

i 

I 


A high  strength  (HS)  layup  consisting  of  90%  0-degree  fibers 
and  10%  90-degree  fibers  within  the  vane  structure.  The 
layup  is  totally  symmetric  with  the  0-degree  fiber  being 
defined  as  a fiber  lying  across  the  length  (8  inch  section) 
of  the  vane.  A 60%  fill  ratio  was  assumed. 

A high  strength  (HS)  layup  consisting  of  10%  0-degree 
fibers  and  90%  90-degree  fibers  within  the  vane  structure. 
(90-degree  rotation  of  test  case  1 layup) 

A high  strength  (HS)  layup  consisting  of  +45-fibers 
throughout  the  entire  vane  thickness.  The  layup  is 
totally  symmetric  and  assumes  a 60%  fill  ratio. 


43 


AFFDL-TR-77-49 


CASE  4l-  Each  of  these  two  cases  is  parallel  to  case  No.'s  1 and  2 ' 

CASE  5) 

respectively.  The  only  exception  is  that  the  high  I 

strength  (HS)  fiber  (Cases  1 and  2)  is  replaced  by  a high  I 

modulus  (HM)  fiber  for  cases  4 and  5.  ] 

CASE  6 - Special  higher  yield  version  of  Case  3. 

CASE  7 - A 21  KSI  fiber  layup  consisting  of  a unit  symmetric  layup 

involving  three  0-degree  layers,  a +45-degree,  a -45-degree  ' 

layer,  and  one  final  90-degree  layer.  Property  data  was 
derived  by  means  based  on  fiber  properties  noted  in 
Reference  6. 

CASE  8 - A 90-degree  rotation  of  the  test  case  7 layup. 

CASE  9 - A compression-molded  vane  of  chopped  filaments;  21  KSI 

fibers  were  assumed  and  orientation  is  random,  resulting 
in  an  isotropic  block  of  vane  material. 

CASE  10  - Baseline  case  on  which  materials  substitution  was  effected. 

Specifically,  this  consists  of  a sintered  metal  matrix 
carbon  graphite  compact  material.  The  material  is 
essentially  isotropic  in  nature  and  can  be  modeled 
several  ways  within  NASTRAN. 

a.  Basepoint  Analysis  (ROVAC,  Case  10)  ' 

Prior  to  conducting  a detailed  materials  substitution  survey,  a 
basepoint  or  baseline  comparison  analysis  was  conducted.  The  basepoint 
analysis  consisted  of  the  ROVAC  graphite  compact  vane  modeled  through  the  | 

use  of  isotropic  NASTRAN  MATI  input  cards.  Several  analysis  checks  and  ; 

debug  runs  were  required  for  each  of  the  associated  load  subcases.  A , 

comparison  of  all  analysis  runs  yielded  good  correlation  and  resulted  in  a j 

final  baseline  run  of  e high  confidence  level.  The  basepoint  analysis  j 

model  consisted  of  an  AVM  two-axle  system  involving  482  gridpoints,  432  i 

quadralateral  elements,  and  67  bar-type  elements.  The  resulting  basepoint  i 

model  is  illustrated  in  Appendix  B Figures  B-4  through  B-6.  A listing  of  j 

the  NASTRAN  input  with  the  exception  of  materials  variation  input  is  j 

contained  in  Appendix  C.  , 

! 

J 
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The  overall  results  of  the  analysis  compared  quite  favorably  with  both 
test  results  and  computations  from  air-cycle  machine  analysis  programs. 
Typically,  bearing  reactions,  deflections,  and  vane  shapes  were  very  plausible 
values  resulting  in  highly  reliable  stress/strain  output  for  the  analysis 
case.  A summary  of  the  results  of  the  MATl  final  run  are  noted  in  the  first 
row  of  Table  5. 

A second  basepoint  run  was  also  made  for  which  a MAT2  input  mode  was 
used.  This  was  done  to  validate  the  MAT2  input  methods  discussed  in  a 
later  section  of  this  report.  The  MAT2  results  correlated  very  well,  as 
can  be  noted  in  Tahle  5. 


TABLE  5 

BASELINE  COMPARISON  SUMMARY 


Lower  Left 
Vane  Deflect 

Vane  Center 
Deflection 

— 

Bearing  Attch 
Deflection 

Top  Center 
Deflection 

Bearing 

Load 

MATl 

BASEPOINT 

.000255 

.00291 

.0000707 

.00298 

971 

MAT  2 

BASEPOINT 

.000256 

.00291 

.0000707 

.00298 

975 

(Deflections  = inches.  Load  = lbs) 


b.  Advanced  Composites  Analysis 

Basically  this  analysis  consisted  of  nine  NASTRAN  runs,  each 
involving  a different  set  of  MAT2  data  which  represented  the  varied 
orthotropic  parameters  for  analysis  cases  1 through  9.  The  first  step  of 
the  analysis  was  to  devise  and  tabulate  the  basic  fiber,  lamina,  and 
laminate  properties  associated  with  the  various  analysis  cases.  This  was 
accomplished  by  extracting  basic  filament  or  specimen  data  (References  1, 

2,  and  4)  and  applying  this  data  to  the  SQ5  program  (Reference  3)  to  obtain 
basic  laminate  properties.  A summary  of  the  extracted  data  and  resulting 
SQ5  properties  is  contained  in  Table  6.  From  this  data  the  basic  bulk  data 
G values  required  by  NASTRAN  can  be  derived  for  MAT2  input.  A complete 
description  of  how  the  MAT2  values  G^j^,  G^j^,  ^22^  ^22'  ^33’ 

derived  and  run  is  included  in  the  following  discussion. 
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of: 


(1)  MAT2  (NASTRAN  Input) 

The  NASTRAN  bulk  data  input  for  the  anisotropic 


The  values  associated  with  the  above  matrix  are 
through  the  Reference  2 equation: 


where: 


n - E (l-y) 

^11  " (l+p)0-2TiT 


Q 


12 


Q.,  = Q 


. P E 


13  ' ^23  ■ {1+vi)(1-2m) 


^44  " ‘"xy 
^55  " Sz 
^66  " ^xz 


and 


'12  ‘■22 


'21  ^11 


case  consists 


(1) 


related 


(2a) 

(2b) 

(2c) 

(2d) 

(2e) 

(3) 
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Q 


The  actual 
values  and 


bulk  data  NASTRAN  conversion  consists  of  Equation  2 

= G . Further  restricting  the  vane  to  an  in-plane 
xy 


G 


13 


= 0, 


from  which  Equation  1 becomes; 

^11  ^21^11 

I-U12U21  ^'^^12^21 

^21^11  ^22 

I-U12U21  ^’1^12^21 

0 0 

(a)  Analysis  Case  1 (MAT2  Input) 

For  the  case  of  [Og*  90]^  and  from  Table  6,  which 
was  prepared  via  composites  program  SQ5  (Reference  3),  the  following 
values  are  extracted:  E^-j  = 19.0(10^),  E22  = 3.6(10^),  G^^  = .65(10^), 
and  ” -0378. 


From  Equations  1 and  4: 

G,^  = (10^)19. 0/[l-(.20)(. 0378)]  = 19.44(10^) 
G^2  = (10^)(.0378)(19. 0)7.99244  = .724(10®) 
G22  = 3.6  (10®)/. 99244  = 3.63(10®) 

*^13  " ^23 
G33  = .65(10®) 
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For  the  ROVAC  AVM  input  the  appropriate  MAT2  input  for  Areas  3 and  4 is 
as  follows: 


ITEM  NO. 

1 

2 

3 

4 

5 

6 

1 

7 

8 

9 

COL.  LOC. 

1-8 

17-24 

25-32 

33-40 

41-48 

49-56 

57-64 

65-72 

VARIABLE 

MAT2 

MID 

Gll 

G12 

G13 

G22 

G23 

G33 

RH0 

Area  3 Inout 

MAT2 

1 

3 

19.14+6 

.724+6 

- 

3.63+6 

- 

.65+6 

.056 

Area  4 Input 

MAT2 

L _ 1 

II 

19.14+6 

.724+6 

- 

3.63+6 

- 

.65+6 

.056 

Area  2 Input 

2 

24.84+6 

4.47+6 

- 

15.54+6 

i 

4.93+6 

1 .146 

Additional  extrapolations  or  computations  were  utilized  within  the  analysis 
for  material  Areas  1 and  2 (Figure  14).  Area  1,  being  a steel  edge  of  the 
vane  boot,  could  be  simulated  by  means  of  a MATl  bulk  data  card.  Since 
this  was  the  case,  no  MAT2  computations  were  required,  because  the  original 
MATl  property  card  was  retained.  For  Area  2,  some  form  of  anisotropic 
property  has  to  be  modeled  to  represent  the  boot  which  houses  the  composite 
vane.  Typically,  this  was  accomplished  by  assuming  a 60/40  ratio  of 
composite-to-steel  structure  from  which  modulus  values  were  computed  and 
checked  by  means  of  the  SQ5  composites  laminate  program.  Actual  differences 
in  modulus  values  between  a direct  proportioning  method  and  SQ5  were  very 
minimal.  The  resultant  output  is  as  follows: 

= 23.55(10^),  E22  = 14.74(10®),  6^^  = 4.93(10®) 

= .288,  and  ^2]  ' -180  (reference  Equation  3), 
from  which  (using  Equations  1 and  4)  we  obtain: 


. 'll 

23.55 

= 24.84(10®) 

11 

1-^12^21 

.948 

_ '22 

14.74 

= 15.54(10®) 

'22 

l-lil2l^21 

" ‘.948 

'13 

= G23  = 0 
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^12  ■ ^21^11 

or  = yi2E22  " ( -29) (14. 1 78) (10®)  = 4.11(10®) 

®12  = ^12/''‘^'12^21  " 4.47(10®) 

C.33  = . 4.93(10®) 

(b)  Orthogonal  Case  (Case  2) 

Areas  3 and  4 (MAT2  (3) (4))  analyses  were  as  follows: 

Analysis  Case  2 consists  of  a 90-degree  rotation  of  Case  1 
layup  and  was  run  to  permit  the  extrapolation  of  findings  between  the  two 
orthogonal  extremes.  For  this  case  the  Table  6 (Col.  A)  directional 
properties  were  reversed  to  obtain  the  following:  E.|.|  = 3.6(10®), 

E22  = 19.0(10®),  Gyy  = .65(10®),  p.,2  = .0378,  and  = -20. 
from  which: 

6^^  = 3.63(10®) 

G22  = 19.144(10®) 

G^2  = (-20)  (3.6)(10®)/.99244  = .725 

®13  " ®23  ^ ® 

G33  = .65(10®) 

Area  2 (MAT2  (2))  analyses  proceeded  as  follows: 

By  simply  reversing  the  Case  1 modulus  values,  we  have: 

E^^  = 14.74(10®),  E22  = 23.55(10®) 

Gj^y  = 4.92(10®) 

P21  = 0.2881  (from  SQ5  reverse) 
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From  this  then; 

= 15.54(10®)  and  0^3  = 24.84(10®) 

0^2  = (0.288)(14.74)  = 4.245(10®) 

^12  ~ ^1 2^^ "'^l 2'^21  ” 4.47(10  ) 

*^13  ""  S3  ^ ° 

S3  ■ Sy  ■ ‘'■«3(loS 

(2)  MAT2  Input  Sunmary  (Cases  1 and  2) 

All  of  the  G values  generated  for  Cases  1 and  2 as  input 
into  the  NASTRAN  bulk  data  set  is  summarized  in  Table  7. 

TABLE  7 

"MAT2"  SUMMARY  CASES  1 AND  2 


VARIABLE 

^ MAT2 

MID 

Si 

S2 

*"13 

G22 

^23 

*"33 

RH0 

& COLUMN 

1-8 

9-16 

17-24 

25-32 

33-40 

41-48 

49-56 

57-64 

65-72 

ANALYSIS 
CASE  NO.l 

MAT2 

2 

24.84+6 

4.47+6 

- 

15.54+6 

- 

4.93+6 

.146 

MAT2 

3 

19.14+6 

.72+6 

- 

3.63+6 

- 

.65+6 

.056 

MAT2 

■■ 

19.14+6 

.72+6 

- 

3.63+6 

- 

.65+6 

.056 

ANALYSIS 
CASE  NO. 2 

MAT2 

15.54+6 

4 +6 

24,84+6 

■ 

4.93+6 

.146 

MAT2 

3 

3.63+6 

.72+6 

1 

19.14+6 

■1 

.65+6 

.056 

MAT2 

■ 

3.63+6 

.72+6 

19,14+6 

■ 

.65+6 

.056 

L 
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(3)  MAT2  Input  Summary  (Cases  1 through  10) 

Following  the  above  procedure  a complete  derivation  of 
all  MAT2  input  requirements  was  then  made.  The  first  step  of  this 
summary  was  to  run  all  of  the  SQ5  runs  required  to  derive  the  laminate 
properties  for  the  40/60  steel -to-composite  section  of  the  vane.  Table  8 
summarizes  all  of  the  values  extracted  from  each  of  the  SQ5  runs. 


The  next  step  was  to  directionally  correct  the  Table  8 values  by 
means  of  Equations  5a  through  5g  into  a set  of  G values  for  each  laminate 
noted  in  Table  8.  The  results  of  this  conversion  are  noted  in  Table  9. 

The  resulting  Table  9 values  were  expanded  to  the  orthogonal  sets  and 
all  of  the  required  Area  2 input  values  were  derived  for  Cases  1 through  8, 
i as  noted  in  Table  10.  Also  contained  within  Tables  9 and  10  are  the 

i associated  values  for  analysis  Cases  9 and  10.  Input  for  these  cases 

j will  be  discussed  in  the  next  section, 

i 

I 

, TABLE  8 

I 

i SQ5  PROPERTIES  FOR  40/60  STEEL  COMPOSITE  BOOT  AREA,  USING  5 DIFFERENT 


COMPOSITE  INSERTS 


HS-FIBER 
[Og  , 90]^ 

HM-FIBER 
(Og  ,90,)^ 

HS-FIBER 

(±45)^ 

SPEC. FIBER 
( !45) 

21  KSl 

AS-FIBER 
(O3.  +45,  90)^ 

Ex 

23.55 

25.73 

14.10 

14.31 

19.44 

'y 

14.74 

15.03 

14.10 

14.31 

15.88 

^xy 

4.93 

4.93 

7.85 

8.42 

5.90 

.288 

.289 

.429 

.450 

.316 

mm 

.180 

.169 

.429 

.450 

.258 

l-M^My 

.948 

.951 

.816 

.798 

.918 

(*Derived  via  Equation  3) 
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The  next  step  in  the  MAT2  summarization  is  the  definition  of  the 
G matrix  for  vane  areas  3 and  4.  This  was  done  through  the  use  of  the 
base  Equations  1 and  4.  From  these  equations  and  data  established  in  Table 
6 the  summarized  results  noted  in  Table  11  for  areas  3 and  4 were  obtained. 


TABLE  9 

"MAT2"  INPUT  G VALUE  AND  ASSOCIATED  EQUATIONS  (X  10®) 


N^YUP 

proKv 

HS-FIBER 

[09,90]s 

HM-FIBER 

[09.90]s 

HS-FIBER 
[ !45]s 

SPEC-FIBER 

[+45]s 

21  KSI 

AS-FIBER 
[O3,  +45,  90]s 

®11 

24.84 

27.05 

17.28 

17.93 

21.18 

Ha 

15.54 

15.80 

17.28 

17.93 

17.30 

mm 

4.47 

4.57 

7.41 

8.07 

5.46 

mm 

0 

0 

0 

0 

0 

mm 

0 

0 

0 

0 

0 

aa 

4.93 

4.93 

7.85 

8.42 

5.90 

Gn  = 

11  ^"^12^^21 

(From  Equation  4) 

(5a) 

^x®y  " ^y^x  ^12®22"  ^21®11 

(5b) 

u h2  ®22  ^x 

(From  Equation  3) 

(5c) 

(5d) 

G , . 2iill  . Mx- 

u l-PxPy  ’"Vy 

(5e) 

®13  " ®23  ^ ° 

(5f) 

®33  = ®xy= 

(5g) 
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(a)  Analysis  Cases  9 and  10 

As  noted  in  Tables  10  and  II,  values  for  two  additional 
cases  were  tabulated.  These  involved  the  chopped-graphite  case  (Case  9) 
and  MAT2  basepoint  run  (Case  10).  Each  of  these  cases  could  have  been 
simplified  to  MATl  NASTRAN  input,  but  were  MAT2  input.  Case  9 was  input 
this  way  for  program  input/output  consistency.  Case  10,  however,  was 
included  for  correlation  to  a NASTRAN  run,  using  MATl  data  cards.  For 
this  correlation  the  MAT2  values  (G^^  G^2>  ^22'  ^23’  ^33^ 

generated  through  the  equations  which  effectively  degenerated  the  MAT2 
data  to  MATl  data  within  NASTRAN.  Correlation  was  excellent  (less  than  1% 
error)  between  the  two  runs,  indicating  a high  level  of  confidence  in  the 
MAT2  input.  The  comparative  results  for  the  Case  10  correlation  were 
noted  back  in  Table  5.  New-data  Cases  9 and  10  input  values  for  Tables  10 
and  11  were  derived  through  Equations  1 and  4 and  the  data  contained  in 
Table  6,  whereby  orthotropic  properties  were  degenerated  by  setting 
Ell  ” ^22'  shear  was  computed  simply  by;  G^^  = E/2(l+p)  with 

U = P^2  " ‘^21- 

(b)  Analysis  MAT2  Input  Summary  For  Cases  1 Through  10 

The  final  step  in  the  MAT2  summarization  was  to  combine 
all  of  the  data  from  previous  tables  into  a summary  table  for  each  of  the 
10  analysis  cases.  This  combination  is  noted  in  Table  12  where  all  of  the 
data  is  ordered  and  set  up  as  to  the  actual  input  format  required  by 
NASTRAN. 

2.  PART  II  (AXLE  AND  BOOT  ELIMINATION) 

The  second  part  of  the  analysis  consisted  of  eliminating  both 
vane  axles  and  the  end  boots  (Figure  4).  These  design  changes  (denoted 
as  FDL-1  and  FDL-2)  were  effected  quite  simply  through  property  changes 
to  those  elements  representing  the  axle  or  boot  portions  of  the  model. 
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VO  VO  ^ VO  VO  VO  VO  voko  VO  VO 
^tovf)  ^lOvowLnvnwiovn 

oo  ooL-oo^oo 


) VO  VO  VO  VO  VO  VO  VO  VO  • 
•mvowvnml^vnvr)* 
•ook-ook-ooi 


>V0V0^V0V0[CSIV0V0 

-vnvn^mvovnvovo 

'OOr-OOb— oo 


U3tnvO  VOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVOVO' 
fovOlO  rovovnl^)<\JCvl^ovoLf>fovov^>CMvovo0^^r*■.or*^^*«.^OOor*^^^ 
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a.  Axle  Removal  (FDL-1) 

Uithin  the  advanced  composite  (Part  I)  model  the  axles  are 
represented  by  7th  and  8th  (Areas  7 and  8,  Figure  14)  MAT  items  in  the 
MAT  data  set.  For  these  items,  the  modules  and  density  properties  were 
assigned  insignificant  values  of  3.0  and  .0001,  respectively.  This,  in 
effect,  yields  a set  of  axles  with  no  weight  nor  stiffness,  thus  eliminating 
them  from  the  model . 

b.  Boot  Removal  (FDL-2) 

'•iithin  the  Part  I model,  the  boots  are  represented  by  the  1st 
and  2nd  (Areas  1 and  2,  Figure  14)  MAT  items  in  the  MAT  data  set.  For 
boot  removal,  the  properties  for  these  areas  were  simply  equated  to  those 
in  areas  3 and  4,  which  resulted  in  a uniform  vane  of  one  material  and 
no  boots. 

3.  PART  III  (VANE  THICKNESS  VARIATIONS)  FDL-2-.xt 

The  third  portion  of  this  study  consisted  of  changing  the  thickness 
of  the  two  most  promising  designs  resulting  from  the  parts  II  and  HI 
studies.  Here  again,  the  model  change  was  quite  simple  because  the 
"PQUAD2’'  data  sets  controlling  the  thickness  valves  in  Areas  1 through  4 
(Figure  14)  were  all  that  needed  revision.  The  specific  analysis  consisted 
of  reducing  the  thickness  to  70%,  60%,  50%,  and  40%  of  the  basepoint 
thickness  values  carried  in  Parts  I and  II. 

4.  PARTS  I,  II,  AND  III  SUMMARY 

A complete  summary  of  all  revisions  to  the  NASTRAN  program  required 
for  the  analysis  is  noted  in  Table  13.  The  data  contained  in  this  table 
was  generated  from  the  actual  NASTRAN  changes  for  each  of  the  22  computer 
runs  comprising  the  analysis  conducted  in  this  report.  Typically,  the 
data  supplied  in  Table  13  were  the  only  changes  made  to  the  NASTRAN 
program  noted  in  Appendix  C. 
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TABLE  13 


ACTUAL  NASTRAN  INPUT  CHANGES 


a.  Analysis  Cases  1 Through  10 


TITLErK.o. 

SrMMARTZ,  ANALTSIS  il 

( 

0 

nrs.  .AYER5=9,  90 

OEG.  LATE 

RS  = 1) 

HS 

1AT2  2 

?4.P4»6 

4.47*5 

15.54*6 

4. 93*6 

. 146 

HAT?  3 

iq.i4*6 

.7?  *5 

3.63*6 

. 65*6 

.056 

HAT?  «. 

IP. 1A»6 

.7?*5 

3.63*6 

. 65  *3 

.055 

TITLF  = »'.‘>.' 

SCHHAPTZ,  ANALYSIS  #? 

r 

0 

OES.  LAY5RS=:1,  90 

DEG.  LAYERS=9» 

HS 

HAT?  2 

4.47*6 

24.84*6 

4. 93*6 

.145 

HAT?  3 

3.63*6 

.33*6 

19. 14*6 

.65*6 

.056 

HAT?  It 

3.63*6 

. 73*6 

19.14*6 

. 65*6 

.056 

TITLf^K.F. 

'JCHHAPT?,  ANALYSIS  *3 

( 

ALTERNtTINC  *45,-45 

DEG.  LAYERS  » 

HS 

HftT2  ? 

17.?A*6 

7.41*6 

17.28*6 

7.85*6 

. 146 

HAT?  3 

6.E3*6 

5.23*6 

6.53*6 

5.52*6 

• 0^3 

HAT?  I* 

6»  5 3 6 

5.23*6 

6.53*6 

5. 52*6 

.055 

TJTLE=tf.o. 

SCHWARTZ,  ANALYSIS  #A 

( 

0 

PEG.  LAYERS=9,  90 

DEG.  LAY; 

RS  = 1) 

HH 

HAT?  ? 

37.C5*6 

4.57*6 

15.8C*  6 

4.93*6 

. 146 

HAT?  3 

?2.8A*6 

1.08*5 

4.15*6 

. 65*6 

.056 

HAT  ? 1. 

??.  9A*6 

1.08*6 

4.15*6 

. 65*6 

.055 

TITL  FrIf.P.' 

FCHHAPTZ,  ANALYSIS  #5 

( 

0 

D'^3.  LAYERS  = 1,  90 

DEG.  LAY; 

RS=9) 

HH 

HAT?  2 

15.70*6 

4.57*6 

27.05*6 

4.93*6 

.145 

HAT?  3 

4.15*6 

1.08*6 

22.88*6 

• 65  **6 

. 055 

HAT?  It 

4.  15*6 

1.08*6 

27.88*6 

. 65*6 

.056 

TITLt  = 'C.P, 

SCHWART7,  ANALYSIS  #6 

( 

ALTERNSTINS  *45,-45 

DEG.  LAY; 

RS  ) 

HH 

HAT?  ? 

17.93*6 

8.07*6 

17,93*6 

8. 42*6 

. 146 

HAT?  3 

3.  45*6 

2.86*6 

3.45*6 

6. 46f 6 

.056 

HAT?  A 

3.45*6 

2.86*6 

3.45*5 

6. 46*6 

. 056 

TITLE^K.P. 

SCHWARTZ,  ANALYSIS  #7 

( 

0 

PEG. 73, *45, -45,90 

DEG.  LAYERS  ) 

AS 

HAT?  ? 

?1.1A*6 

5.46*6 

1 7. 30*6 

5.  90*6 

.146 

HAT?  3 

l’.PA*6 

1.97*6 

6.53*6 

2. 27*6 

• 0 56 

HAT?  It 

1 3.PA*6 

1.97*6 

6.53*6 

2.27*6 

.056 

title=ii.?. 

SCHWARTZ,  AI'ALYSIS  #8 

193 

PEC. =3, *45, -45,  0 

DFG.  LAYERS  ) 

AS 

HAT?  ? 

17.30*6 

5 • 46 ♦6 

21.18*6 

5. 90*6 

. 146 

HAT?  3 

5.53*6 

1.98*6 

12.98*6 

2. ?7*6 

.055 

HAT?  4 

5.53*6 

1.98*6 

12. 98*6 

2.27*6 

.056 

TITLE=K.P. 

SCHWARTZ,  analysts  *9 

{ 

CHOPPED  GRAPHITE  HATERIAD 

HAT  ? ? 

20.05*6 

5.99*6 

20.05*6 

7. 74*6 

.146 

hat?  3 

14.82*6 

4.42*6 

14.B’*t 

5.  20*6 

. 055 

HAT?  A 

14. 8?*6 

4.4?*b 

14.82*6 

5.20*6 

. 055 

TITL  E^H.o. 

SCHWARTZ,  analysis  *10 

{ 

GRAPHITE  COHPACT  HATERIAL) 

“AT?  ? 

15.1?*6 

3.67*6 

15. 12*6 

5. 60*6 

.152 

HAT?  3 

3.60  *6 

1.06*6 

3. 80*6 

1. 37*6 

.066 

HAT?  A 

3.  80*6 

1.06+6 

3. 80*6 

1. 37*6 

. 066 
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TABLE  13  (CONI' D) 

ACTUAL  NASTRAN  INPUT  CHANGES 
b.  Analysis  Cases  7A  Through  70 


TITLF^K. 

CHMAPT7,  NO  axles  #7 

( 

J npj,  = 3f*45f-45|'?0 

PEG.  LAfiRS  ) 

A> 

MAT? 

? 

?1. 1"*6 

5.96*6 

1 7.33*6 

5«  Q 0 *-6  • 146 

HAT? 

1?,3A*6 

1.97*6 

6.53*6 

2.27*6  .056 

MAT? 

A 

1?.'>8*B 

1.97*6 

6.53*6 

2.27*6  .056 

HATl 

7 3.0EC 

.32c:  .ccri 

“ATI 

8 3. CEO 

.3280  . 03  31 

TITLE^tr. 

P.SCHHA^T-r,  NO  BOOTS  #73 

( 

0 0E3. =3,*95,-95,90 

OEG.  LATEPS  3 

AS 

MAT2 

1 

l?.a8K> 

1.07*6 

6.53*6 

2.27*6  .056 

mat? 

? 

1?.98*6 

1.97*6 

6.53*6 

2.  27*6  .056 

HAT? 

T 

12.98*6 

1.97*6 

6.53*6 

2.27*6  .055 

MAT? 

U 

1?.°8»6 

1.97*6 

6.53*6 

2.  27*6  . 055 

MATl 

7 ■'.CEO 

.i2:c  .cci 

“ATI 

8 3.0EO 

. 3200  . 0031 

TITLf  =f. 

».SCHHART7,  ?Z7.  size  #70 

( 

0 OE*.. =3, *95, -95, 90 

DEG.  LATERS  > 

A3 

PqUAP? 

1 1 

.2620 

PQUAC? 

7 ? 

.2620 

POUAP? 

3 3 

.2627 

^OUAP? 

4 4 

.2620 

MAT? 

1 

12.98*6 

1.97*6 

6.53*6 

2.27*6  . 055 

HAT? 

? 

12.98*0 

1.97*6 

6.53*6 

2.  27*6  .056 

HAT? 

3 

12.98*6 

1.97*6 

6.63*6 

2.  27*6  . 056 

MA^a 

1?,Q8*-C 

1.97*6 

6.53*6 

2.27*6  . 055 

HAT! 

7 3.0E0 

. 3203  .0031 

HATl 

3 3.0EC 

.3203  .0381 

TITLE^tf. 

'’.SCHWf^r/,  60,;  SI7E  ,73 

( 

0 nEj»=’3ffr4';,-45f5!) 

OEG.  LATEPS  t 

AS 

POUAO? 

1 1 

.2290 

PqUAT? 

2 2 

.2290 

pqijAc? 

3 3 

.2290 

P^rtr  C ^ 

4 

.229 

MAT? 

1 

12.98*6 

1.97*6 

6.53*6 

2.  27*6  . 056 

HAT2 

9 

l’.98*6 

1.07*6 

6.53*6 

2.  27*6  . 055 

MAT? 

7 

12.98*6 

1.97*6 

6.57*6 

2.  27*6  .056 

HAT? 

4 

12.98*6 

1.97*6 

6.53*6 

2.  27*6  . 056 

“ATI 

7 3.OE0 

.3200  .O'Cl 

«AT1 

8 3.0EO 

. 3203  . 00"! 

c. 

Analysis  Cases  7E  and  7F 

TITLP  = «’. 

O.SCHWA3T7,  SC'  SIZE  #7E 

0 PEG. =3, *95, -95, 90 

DEG.  LATERS  » 

AS 

pouA  n? 

1 1 

.1870 

POUAf? 

? ? 

.1870 

POUAP? 

3 3 

.187.1 

PQUAP? 

4 4 

.18’7 

HAT  ? 

1 

1’.98*6 

1.97*6 

6,53*6 

2.27*6  . 055 

HAT? 

? 

12.98*6 

1,97*6 

6.53*6 

2.  27*6  . 056 

MAT? 

T 

1 2.98*6 

1.97*6 

6.53*6 

2.27*6  .056 

MAT  ? 

U 

12.98*6 

1.97*6 

6.5-'*6 

2.  27*6  . 056 

HATl 

’ 3. CEO 

.3200  .00"! 

HATl 

• 3.0E0 

. 3 203  .OCl 

TITLE^t 

.O. 

SCHhAoTZ,  903 

SI3E  #7- 

( 

0 nSG. =3,*95,-95,93 

DEG.  LATERS  ) 

AS 

POUAP? 

1 1 

. 1500 

POUAP? 

2 2 

. 150C 

POUAP? 

3 3 

.15oI 

PQUA  n? 

4 4 

. 1=00 

MAT  ? 

1 

12.98*6 

1.97*6 

6.53*6 

2.  27*6  . 055 

HAT? 

? 

12.08*6 

1.97*6 

6.53*6 

2.  27*6  . 055 

HAT? 

7 

12.98*6 

1.97*6 

6.57*6 

2.  27*6  .055 

MAT  ? 

4 

12.08*6 

1.97*6 

e.5'*6 

2.  27*6  . 055 

MATl 

7 3.0EO 

.3330  •C?’’! 

HATl 

8 3.0E0 

. 3 200  .C"’’! 
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TABLE  13  (CONT'D) 


TITLED''. 

P. 

?CHWA?T?,  MC  A)(L"S 

CT0PP£0  G94PMITE 

MATERIAL) 

HAT  ? 

? 

Ri.riTfc 

5.99*6 

20.05*6 

7.74*6  .146 

MAT? 

1I*.S2A6 

4.42*6 

14.82*6 

5.20*6  .056 

MAT  ? 

<A 

lA. ?2»6 

4.42*6 

14.82*6 

5.20*6  .085 

MAT  1 

7 3. "EC 

.3203  .0001 

MATl 

i 3. CEO 

. 3200  . 0001 

TITLE  = K’. 

P. 

5CHMAPT?,  NO  BOOTP 

CMOPPTD  graphit; 

MATERIAL) 

MAT? 

1 

1A.A2T9 

4.42*6 

14.  82*6 

5.  20*6  .055 

MAT? 

? 

1A.8?"& 

4.42*6 

14.  82*  6 

5.20*6  .056 

MAT? 

T 

1A.8?"*; 

4.42*6 

14.82*6 

5.  20*6  . 056 

MAT  2 

4 

14.82TP 

4.4?*6 

14.82*6 

5.20*6  .056 

'lATl 

7 3. CEO 

.32Cj  .DTOl 

MATl 

« 3. CEO 

.32C0  . 000  1 

TITLF-if. 

p. 

FCHHARTZ,  707 

SIZE  #92 

CTOP^En  GRAPHIIi 

MATERIAL) 

oOOAn? 

1 1 

.2623 

POUAPB 

2 2 

.2620 

PQUA  n3 

T 3 

.2623 

POUAP? 

4 4 

.2620 

HAT  2 

1 

14.C2T6 

4.42*6 

14.82*6 

5.20*6  .056 

MAT? 

? 

lA.A’TE 

4.42*6 

14.82*6 

5.20*6  .055 

MAT  ? 

1 

14.  8?*C- 

4.42*6 

14.82*6 

5.  20*6  .056 

HAT? 

4 

14.e2T6 

4.42*6 

14.82*6 

5«  ?C  + 6 •056 

MATl 

7 3. CEO 

. 3 203  .COOl 

MATl 

8 3. CEO 

. 3203  .COCl 

TITLE  = »'. 

p. 

3CHMAPTZ,  805: 

site  #90 

CMOPPEP  GRA"MITE 

MATERIAL) 

roUA  D2 

1 1 

.2243 

POUA  n-J 

’ 2 

.2240 

P0M4n2 

T 3 

.2?40 

POUAP’ 

4 4 

.224 

MAT? 

1 

14.82»6 

4.42*6 

14.82*6 

5.  20*6  . 056 

MAT? 

? 

14.82*6 

4.42*6 

1 4. 82*6 

5.  20*6  . 056 

MAT? 

3 

14.82*6 

4.42*6 

14.82*6 

5«  ZZ ^6  *056 

MAT? 

4 

14.82*6 

4.42*6 

14.82*6 

5.20*6  .055 

MAT! 

7 3. CEO 

.32C0  .OO"! 

MATl 

8 3.0EC 

. 3 200  . 0001 

e.  Analysis  Cases 

9E  and  9F 

TITLE:"'. 

®*SCHWART7i  ^IZE 

< 

CTOPPEO  GOAOHITI 

MATERIAL) 

PQUAP? 

1 1 

.1670 

POUAD? 

2 2 

.1870 

poufln? 

T 3 

.1870 

POUAP? 

4 4 

.1870 

MAT? 

1 

14.82*6 

4.42*6 

14.8?*6 

5.  20*6  .055 

HAT3 

? 

14.82*6 

4,42*6 

14.82*6 

5.20*6  .056 

MAT? 

3 

14.87*6 

4.42*6 

14.82*6 

5.  20*6  , 056 

MAT? 

4 

14.8Tf6 

4.42*6 

14. "2*6 

5a  ?0**6  • 05o 

MATl 

7 3. PEC 

.3200  .0001 

MATl 

8 3.CEC 

. 3200  .COOl 

TITLE:"'. 

SCmmAR:?,  407.  SIZE  *9" 

< 

CHOPPEO  GRAPHITE 

material) 

POUAP? 

1 1 

.1503 

POMAO'^ 

2 2 

. 1500 

POUAO? 

3 3 

. 1500 

POUAD? 

4 4 

. I'OO 

MAT? 

1 

14.82*6 

4.42*6 

1 '*.82*6 

5.  20*6  . 056 

MAT? 

? 

14.82*6 

4.42*6 

14.82*6 

5.  20*6  .056 

MAT? 

3 

14.87*6 

4.42*6 

1 4. 82*  6 

5.20*6  .055 

mat  ? 

u 

14.8?*6 

4.42*6 

14. 82*6 

5.20*6  . 055 

MATl 

T 3.060 

,3203  .COCl 

MATl 

8 3. CEO 

. 3 203  .COCl 
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SECTION  V 
VANE  ANALYSIS 

The  overall  analysis  consists  of  a total  of  22  NASTRAN  runs  summarized 
in  Table  14.  The  first  of  these  runs  consists  of  a single  baseline  run 
(Case  10,  Table  14)  representing  the  current  ROVAC  vane.  The  next  series 
of  runs  consists  of  nine  surveys  (Cases  1 through  9,  Table  14)  investi- 
gating several  different  layups  and  fiber  types  with  which  a ROVAC  vane 
could  be  constructed.  Two  of  the  nine  cases  were  selected  as  the  most 
promising  for  further  analysis. 

The  specific  cases  selected  for  further  study  were  Cases  7 and  9, 
primarily  due  to  the  minimal  tip  deflections  experienced  with  these  two 
designs.  The  additional  analysis  consisted  of  first  removing  the  axles 
from  the  vane  (subcase  A)  and  secondly  removing  both  the  axles  and  boots 
from  the  vane  (subcase  B).  These  cases  are  referred  to  as  design 
types  FDL-1  and  FDL-2,  respectively,  and  are  noted  as  test  cases  7A,  76, 
9A,  and  9B  in  Table  14. 

Following  the  axle  and  boot  removal,  subcases  C,  D,  E,  and  F were 
run  on  FDL-2  design  and  consisted  of  thinning  the  FDL-2  design  to  .7, 

.6,  .5,  and  .4,  of  the  original  thickness.  These  cases  are  referred  to 
by  their  thickness  value  (t)  added  to  the  design  type  designation.  The 
actual  analysis  cases  are  noted  as  test  cases  7C  through  7F  and  9E 
through  9F.  An  additional  designation  (CG)  has  also  been  added  to  the 
design  type  to  distinguish  the  chopped-graphite  material  from  the 
layup  type  structures  of  Case  7. 

Within  the  analysis,  the  NASTRAN  case  control  program  (Appendix  C, 

Pg  120)  was  written  to  yield  node  point  displacements,  load  vectors, 
single  point  constraint  forces,  quadrilateral  and  bar  element  forces, 
and  quadrilateral  and  bar  element  stresses.  In  addition  to  this  tabular 
output,  two  plots  consisting  of  single  and  superimposed  deflected  shapes 
were  also  programmed. 
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1.  PART  I ANALYSIS  AND  DISCUSSION 

Actual  data  extracted  from  the  NASTRAN  output  was  limited  to  six 
primary  deflections,  the  ten  highest  stresses  in  the  vane  itself,  and 
the  resulting  bearing  loads  and  the  vane  weight.  The  six  primary 
deflections  of  interest  are  noted  in  Figure  15.  Actual  quantitative 
deflection  values  (Figure  15)  for  Cases  1 through  10  are  included  in 
Table  15,  along  with  vane  weight  and  resulting  bearing  loads.  For  this 
analysis  and  the  following  Part  II  and  III  studies,  deflection,  load, 
and  weight  data  were  extracted  as  follows: 

Deflection:  A = T1  or  X deflection  0 Node  point  1 

B = T2  or  Y deflection  0 Node  point  1 

C = T2  or  Y deflection  O Node  point  28 

D = T2  or  Y deflection  0 Node  point  309 

E = T2  or  Y deflection  0 Node  point  476 

Bearing  Load:  Summation  of  T2  SPC  forces  at  node  points  477,  478, 
479,  and  480. 

Vane  Weight:  2X  mass  noted  in  grid  point  weight  generator. 

Element  and 

Nodal  Locations:  Refer  to  pages  101,  102,  103,  Appendix  B 

The  resultant  stress  data  is  noted  in  Tables  16a  and  16b.  This 
tabulation  notes  the  ten  highest  stresses  occurring  in  the  vane.  All  of 
the  stresses  noted  in  Table  16  occur  in  either  the  vane  boot  region  or 
the  vane  closeout  structure.  In  view  of  this  fact.  Tables  17a  and  17b 
are  also  included,  in  which  is  related  detailed  stress  data  as  restricted 
to  the  composite  portion  of  the  vane  only. 

For  the  most  part,  the  data  presented  in  Tables  15,  16,  and  17  is 
self-explanatory.  Tip  deflection  for  cases  1 through  9 ranged  from 
.93  mils  to  3.74  mils,  as  compared  to  the  basepoint  ROVAC  value  of 
2.91  mils.  It  may  be  well  to  note  that  this  2.91  basepoint  value 
compares  favorably  with  an  estimated  3.0  value  assessed  in  an  AFFDL 
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TABLE  16 

MAXIMUM  PRINCIPAL  STRESS  SUMMARY  REFLECTING  THE  TOP  TEN  STRESS  VALUES 
a.  Analysis  Cases  1 Through  4 and  Baseline 


BASELINE 

ANALYSIS  NO.  1 

ANALYSIS  No.  2 

ANALYSIS  No.  3 1 

ANALYSTS  No.  4 i 

ELEMENT 

STRESS 

ELEMENT 

STRESS 

ELEMENT 

STRESS 

ELEMENT 

STRESS 

ELEMENT 

STRESS 

163 

3035 

163 

2465 

246 

3031 

163 

2861 

163 

2426 

190 

4679 

190 

4158 

190 

3921 

190 

4023 

190 

4093 

191 

3908 

191 

3553 

191 

3254 

191 

3223 

191 

3501 

217 

3887 

215 

. 

3773 

217 

3397 

217 

3115 

217 

3728 

218 

3906 

218 

3670 

218 

3351 

218 

3118 

218 

3626 

244 

4931 

244 

5111 

244 

4402 

244 

3805 

244 

5070 

245 

4838 

245 

4275 

246 

4206 

245 

3937 

245 

4219 

271 

8608 

271 

6241 

271 

7697 

271 

7649 

271 

6086 

272 

7004 

272 

4468 

272 

6212 

272 

6405 

272 

4334 

273 

3160 

164 

2594 

273 

3110 

273 

3240 

164 

2554 

Stress 

in  PSI 

b. 

Cases  5 

Through  9 

ANALYSIS  Ho.  5 

ANALYSIS  No.  6 

ANALYSIS  No.  7 

1 ANALYSIS  No.  8 

1 ANALYSIS  No.  9 1 

ELEMENT 

STRESS 

ELEMENT 

STRESS 

ELEMENT 

STRESS 

IjgggH 

STRESS 

246 

3152 

163 

2975 

163 

2462 

163 

2395 

244 

4572 

190 

3863 

190 

3969 

190 

3917 

190 

3800 

190 

3571 

191 

3208 

191 

3101 

191 

3346 

191 

3208 

191 

3038 

217 

3399 

217 

2821 

217 

3504 

217 

3348 

217 

3240 

218 

3336 

218 

2870 

218 

3492 

218 

3349 

218 

3249 

244 

4458 

244 

3209 

245 

4204 

245 

4147 

245 

4027 

245 

4182 

245 

3772 

246 

2515 

246 

2749 

246 

2668 

271 

7564 

271 

8269 

271 

6436 

271 

6985 

271 

6259 

272 

6039 

272 

7356 

272 

4780 

272 

5435 

272 

4688 

273 

3122 

273 

3670 

1 

164 

2480 

273 

2327 

273 

2617 

Stress  in  PSI 
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TABLE  17 

COMPOSITE  REGION  STRESS  SUMMARY 
a.  Analysis  Cases  1 Through  10 


ELEMENT 

NORMAL 

NORMAL 

SHEAR 

MAOOR 

ANGLE 

MAX. 

X 

Y 

XY 

PRINCIPAL 

a 

SHEAR 

ft  " 

27 

1875 

13.1 

-4.5 

1875 

-.14 

81 

1586 

44.8 

-15.4 

1586 

-.57 

■B 

vr 

227 

1 

450 

97.0 

745 

1039 

-38.34 

■9 

25A 

414 

22.7 

704 

949 

-37.22 

K EiS^I 

281 

321 

-74.2 

-569 

726 

-35.43 

308 

240 

-111.7 

-411 

511 

-33.42 

IB 

* 

27 

2313 

8.95 

5.28 

2313 

-.13 

■ 

81 

1558 

35.4 

14.2 

1558 

-.54 

■■ 

** 

227 

435 

0 

-556 

815 

-34.34 

597 

254 

363 

-67.5 

-507 

700 

-33.49 

552 

281 

178 

-431 

-398 

451 

-34.36 

427 

308 

-115 

-113 

-325 

211 

-45.09 

325 

* 

27 

1684 

10.63 

3.8 

1684 

-.13 

837 

81 

1366 

39.32 

-13.70 

1366 

-.59 

663 

ft* 

227 

367 

-137 

-559 

728 

-32.84 

613 

254 

306 

-352 

-511 

585 

-28.60 

608 

281 

203 

-642 

-415 

373 

-22.23 

592 

308 

124 

-486 

-306 

251 

-22.54 

432 

ft 

27 

1386 

18.88 

-4.9 

1386 

B 

683 

81 

1225 

43.81 

-14.63 

1226 

BB 

591 

ftft 

227 

439 

299 

-730 

1103 

-42.25 

734 

254 

443 

282 

-726 

1093 

-41 .84 

730 

281 

386 

211 

-631 

936 

-41.04 

637 

308 

322 

152 

-494 

738 

-40.12 

502 

27 

2373 

8.82 

-5.52 

2373 

-.14 

!!■ 

81 

1548 

34.78 

-14.45 

1549 

-.55 

■9 

ftft 

227 

439 

9.07 

-553 

817 

-34.35 

593 

254 

369 

-63.8 

-505 

702 

-33.38 

549 

281 

170 

-134 

-397 

443 

-34.55 

425 

308 

-147 

-108 

-330 

203 

-46.67 

331 

*The  two  highest  stress  values. 
♦♦Bearing  support  region  stress. 
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TABLE  17  (Cont'd) 
COMPOSITE  REGION  STRESS  SUMMARY 
b.  Analysis  Cases  5 Through  9 


ELEMENT 


NORMAL 

X 


NORMAL 

Y 


SHEAR 

XY 


MAJOR 

PRINCIPAL 


ANGLE 

a 


MAX. 

SHEAR 


27 

81 


1741 

1402 


10.84 

39.48 


-3.86 

-13.56 


1741 

1402 


-.13 

-.57 


865 

681 


227 

254 

281 

308 


376 

312 

199 

110 


-175 

-412 

-722 

-524 


-556 

-507 

-411 

-304 


721 

572 

356 

231 


-31.81 

-27.22 

-20.87 

-21.89 


621 

623 

617 

438 


27 

81 


227 

254 

281 

308 


1008 

902 


21.54 

40.74 


6.95 

19.8 


1008 

903 


294 

325 

315 

315 


230 

246 

217 

195 


-680 

-705 

-639 

-579 


944 

992 

908 

837 


-.404 

1.31 


-43.63 

-43.40 

-42.81 

-42.06 


493 

431 


681 

706 

641 

582 


27 

81 


2007 

163/ 


10.36 

39.14 


-3.82 

-12.14 


2007 

1637 


-.11 

-.44 


998 

799 


227 

254 

281 

308 


512 

478 

349 

133 


87.4 

9.04 

-66.4 

-63.5 


-674 

-611 

-474 

-324 


1007 

899 

659 

374 


-36.25 

-34.5 

-33.18 

-36.6 


706 

655 

518 

339 


27 

81 


1818 

1529 


11.1 

40.5 


-3.52 

-11.93 


1818 

1529 


-.11 

-.46 


903 

744 


227 

254 

281 

308 


480 

449 

347 

213 


102 

-17 

-146 

-150 


-683 

-621 

-485 

-315 


1001 

880 

645 

396 


-37.28 

-34.71 

-31.54 

-30.03 


709 

664 

545 

364 


27 

81 


1980 

1665 


11.23 

40.50 


-3.54 

-11.36 


1980 

1665 


-.10 

-.40 


984 

812 


227 

254 

281 

308 


579 

549 

415 

182 


239 

139 

40 

20 


-730 

-653 

-501 

-316 


1160 

1029 

762 

427 


-38.44 

-36.30 

-34.73 

-37.77 


750 

685 

534 

326 


♦The  two  highest  stress  values. 
♦♦Bearing  support  region  stress. 
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ROVAC  test  effort.  In  general,  stress  for  all  of  the  design  cases  1 
through  9 were  below  that  of  the  basepoint  case.  Since  the  composite 
designs  of  interest  are  far  superior  to  the  basepoint  design,  from  a 
structural  point  of  view,  no  further  study  was  made  relative  to  the 
resultant  stresses.  As  noted  in  Table  15,  vane  weight  was  reduced 
from  .942  lbs  to .864  lbs.  Due  to  this  reduced  mass,  not  only  were  the 
stresses  reduced  but  the  bearing  load  at  the  axle  support  was  reduced 
from  975  lbs  to  895  lbs. 

Based  primarily  on  the  maximum  tip  deflection  (value  C,  Table  15), 
designs  from  Cases  7 and  9 were  selected  for  further  study.  Specifically 
for  these  cases,  tip  deflections  were  reduced  from  2.91  miles  to  1.20 
and  .93  mils.  One  secondary  point  to  note  is  that  the  movement  toward 
the  compressor  end  plate  (A,  Figure  15)  was  also  reduced  by  over  60%. 

This  could  very  well  contribute  to  a quieter  and  cooler  running  system 
due  to  reduced  end  plate  forces. 

Deflected  shapes  for  all  of  the  analysis  cases  (Parts  I,  II,  and  III) 
are  included  in  Appendix  D.  These  plots  include  the  deflected  shape 
superimposed  on  the  original  static  shape.  It  should  be  noted  that 
Appendix  D figures  cannot  be  compared  to  each  other  because  NASTRAN 
plotting  equates  the  maximum  deflection  to  a specific  value. 

2.  PART  II  ANALYSIS  AND  DISCUSSION 


For  the  second  analysis,  the  primary  items  of  interest  are  the  same 
as  for  the  previous  case.  The  overall  bearing  load,  vane  weight,  and 
deflection  summary  data  are  noted  in  Table  18.  Changing  regions  of 
maximum  stress  can  be  noted  in  Table  19  and  quantitative  stress  data 
is  included  in  Table  20.  Again,  element  and  nodal  locations  can  be 
obtained  by  referring  to  Appendix  B (Pages  101,  102,  103).  All 
deflection,  loads,  and  weight  parameters  are  the  same  as  in  Part  I. 


The  data  has  again  been  tabulated  in  a very  self-explanatory  manner. 
Tip  deflections  tended  to  increase  about  14%  for  the  FDL-1  (no  axle) 
design.  For  the  FDL-2  (no  axle,  no  boot)  design,  additional  tip  deflection 


L 


70 


AFFDL-TR-77-49 


TABLE  18 

LOAD/DEFLECTION  CHANGES  DUE  TO  VANE  REDESIGN 


FIG. 

DEFLECTION 

VALUE 

ANALYSIS 

No.  7 

anaLVsTS' 

No.  7A 
AXLES 
REMOVED 

ANALYSIS 

No.  7B 

AXLES  & BOOTS 

REMOVED 

— 

ANALYSIS 
No.  9 

j 

ANALYSIS 
No.  9A 
AXLES 
REMOVED 

SHalysis 

No.  9B 

AXLES  & BOOTS 
REMOVED 

A 

.00041 

.00043 

.0004 

.00037 

.00039 

.00037 

B 

.0002 

.00035 

.00043 

.00019 

.00034 

.00031 

C 

.00120 

. COl 40 

.00145 

.00093 

.00111 

.00106 

D 

.00005 

.0002 

.00018 

. 000029 

.00022 

.00017 

E 

.00123 

.0014 

.0015 

1 

.00094 



.00112 

.00106 

C-B 

.001 

.00105 

.00102  j .00074 

4-  - - . - 4 

.00077 

.00075 

BEARING 

LOAD  (lbs) 

894 

801 

626 

894 

801 

1 

626 

VANE 

WEIGHT(lbs) 

1 

1 

.864 

.772 

.608 

.864 

.772 

.608 
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TABLE  19 

ELEMENT  NUMBERS  DENOTING  THE  TEN  HIGHEST  MAXIMUM  PRINCIPAL  STRESS  VALUES 
FOR  EACH  ANALYSIS  CASE 
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TABLE  20 

MAXIMUM  PRINCIPAL  STRESS  SUMMARY  FOR  VANE  REDESIGN 
(INCLUDES  ALL  TABLE  ELEMENTS) 


i - ■ ' u 

'“v^NALYSIS 

FI  FMm^ 

No.  7 

No.  7A 

No.  7B 

No.  9 

No.  9A 

No.  9B  1 

1 

24 

1900 

2000 

1866 

1879 

1973 

1851  ; 

25 

1954 

1930 

2023 

1899  i 

26 

1989 

2088 

1957 

1963 

2055 

1931 

27 

2007 

2106 

1975 

1980 

2071 

1947 

163 

2462 

2540 

1372 

2568 

2331 

1307 

164 

2480 

2580 

1320 

2237 

2326 

1307 

190 

3917 

4165 

2659 

3571 

3793 

3034 

191 

3346 

3566 

1946 

3038 

3250 

3962 

217 

1 

3504 

4182 

1978 

3240 

3904 

2262 

218 

3492 

3886 

1893 

3249 

3677 

2198 

219 

2183 

2464 

1753 

2244 

2564 

2022 

244 

4810 

6175 

2425 

4572 

5965 

3277 

245 

4204 

4081 

1994 

4027 

4010 

2462 

246 

2048 

2668 

2597 

2048  1 

n 

271 

6436 

2122 

1455 

6259 

2085 

1582 

272 

^780 

1739 

1360 

4688 

1784 

1460 

273 

2404 

1424 

1260 

2617 

1493 

1327 

D 

Stress  Values  in  PSI 
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changes  were  insignificant.  No  appreciable  change  was  noted  over  the 
already  improved  lateral  movement  toward  the  compressor  end  plates.  In 
general,  no  appreciable  changes  occurred  to  overall  deflections  in 
either  the  FDL-1  and  FDL-2  cases. 

In  removing  the  vane  axles  and  boots,  vane  weights  were  reduced 
further  to  .772  and  .608  lbs,  thereby  substantially  improving  the  bearing 
load  properties.  In  addition  to  the  bearing  load  improvements,  maximum 
principal  stresses  were  also  reduced  for  the  most  part  by  almost  50%  for 
the  FDL-2  case.  For  the  FDL-1  case,  a 5 to  10%  increase  in  stress  values 
was  noted. 

In  summary,  it  can  be  said  that  removing  the  axles,  or  axles  and 
boots,  does  little  to  degrade  the  basic  advanced  composite  approach. 
Improvements,  however,  are  substantial  in  the  area  of  load  and  stress 
reduction,  particularly  for  the  FDL-2  case.  Based  on  these  factors,  the 
FDL-2  design  was  selected  for  further  analysis  in  the  Part  111  section 
of  this  study. 

3.  PART  III  ANALYSIS 

The  Part  III  analysis  consisted  of  defining  the  effects  of  varying 
thickness  on  the  FDL-2  design  selected  in  the  previous  section.  Overall 
deflection  results  for  this  analysis  are  included  in  Table  21  (FDL-2)  and 
Table  22  (FDL-2(CG)).  Stress  data  is  tabulated  in  Tables  23  and  24  for 
the  same  respective  designs.  The  primary  item  to  note  from  this  analysis 
is  that  although  the  amount  of  structure  is  reduced  substantially,  the 
effects  upon  vane  stresses  are  only  minor.  In  addition,  with  the 
decreased  structure,  the  resultant  deflections  of  interest  were  actually 
improved  upon  slightly.  The  prime  improvement,  however,  was  in  the  area 
of  vane-bearing  loads,  where  values  ranged  from  626  lbs  down  to  328  lbs 
for  a vane  of  40%  of  the  original  thickness. 
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TABLE  21 

LOAD/DEFLECTION  CHANGES  DUE  TO  VANE  RESIZING  (MATERIAL  = 0 Deg.  = 3,  +45, 

-45,  90  Deg.)  AS  Type 


FIG. 

DEFLECTION 

ANAL/7B 

ANAL/7C 

ANAL/70 

ANAL/7E 

ANAL/7F  1 

VALUE 

t=100% 

t=70% 

t=60% 

t=50% 

t=40%  1 

A 

.00040 

.00040 

.00040 

.00040 

.00040 

B 

.00043 

.00038 

.00036 

.00035 

.00033 

C 

.00145 

.00139 

.00137 

.00134 

.00132  1 

0 

.00018 

.000145 

.000131 

.000117 

.000103 

E 

.0015 

.0014 

.00138 

.00136 

.00133 

C-B 

.00102 

.00102 

.00101 

.00099 

.00099 

BEARING 

LOAD  (lbs) 

626 

477 

428 

379 

328 

VANE 

.608 

.465 

.417 

.370 

.323 

WEIGHT  (lbs) 

i y 

TABLE  22 

LOAD/DEFLECTION  CHANGES  DUE  TO  VANE  RESIZING  (MATERIAL  = CHOPPED  GRAPHITE) 


DEFLECTION 

VALUE 

r 

ANAL/9B 

t=100% 

ANAL/9C 

t=70% 

ANAL/9D 

t=60% 

ANAL/9E 

t=50* 

l~"  ~ 

ANAL/9F 

t=40% 

A 

.00037 

— 

.00037 

.00036 

.00036 

.00036 

B 

.00031 

.00026 

. 00024 

.00023 

.00021 

C 

.00106 

.00100 

.00097 

.00095 

.00093 

D 

.00017 

.00017 

.00012 

.00011 

.00010 

E 

.00106 

.00100 

.00099 

.00096 

. 00094 

C-B 

.00075 

.00074 

.00073 

.00072 

.00072 

BEARING 

LOAD  (lbs) 

626 

477 

428 

379 

328 

VANE 

WEIGHT  (lbs) 

.608 

.465 

.417 

.370 

.323 

c r 
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TABLE  23 

MAXIMUM  PRINCIPAL  STRESS  SUMMARY  FOR  CHANGES  DUE  TO  VANE  RESIZING 
(TEN  HIGHEST  VALUES)  (MATERIAL  = 0 DEG.  = 3,  +45,  -45,  90  DEG.) 


■^ALYSIS 

.emenYs^ 

No. 

7B 

Mo.  7C 

No.  7D 

No.  7E 

1 

No.  7F 

24 

18 

166 

1856 

1850 

1844 

1835 

25 

19 

121 

1910 

1905 

1898 

1890 

26 

19 

157 

1946 

1940 

1933 

1924 

27 

19 

175 

1963 

1957 

1951 

1942 

190 

26 

.59 

2780 

2826 

2874 

! 2922 

191 

19 

146 

1991 

2007 

2025 

2042 

217 

19 

178 

1940 

1923 

1901 

1876 

218 

18 

193 

1846 

1825 

1800 

1768 

244 

2i 

125 

2213 

2130 

2040 

1938 

245 

19 

194 

1870 

1819 

1 

1762 

' 1697 

1 

TABLE  24 


MAXIMUM  PRINCIPAL  STRESS  SUMMARY  FOR  CHANGES  DUE  TO  VANE  RESIZING 
(TEN  HIGHEST  VALUES)  (MATERIAL  = CHOPPED  GRAPHITE) 
'•^ALYsis  I I ■ 

No.  9B  No.  9C  No.  9D  No.  9E  No.  9F 

ELEMENr\ 
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The  primary  factor  to  be  derived  from  this  analysis  is  that  for  this 
i particular  design  an  optimum  design  would  be  the  thinnest  vane  possible. 

' Obviously  other  factors  come  into  play  such  as  manufacturing  limitations 

; and  out-of-plane  considerations.  Typically,  these  considerations  would 

, include  vane  bending  along  with  the  complex  aspects  of  modeling  the 

I structure  with  3-dimensional  isoparametric  elements. 

1 

I 

I 

1 

I 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

The  results  of  the  programs  and  analysis  conducted  within  this 
program  has  led  to  the  following  observations  or  conclusions: 

1.  For  the  specific  design  considered  within  the  analysis  (Figure  4), 
the  existing  carbon  compact  type  material  can  be  improved  upon  substantially 
through  the  use  of  advanced  composites.  These  improvements  can  further 
lead  to  more  permissive  design  variation  for  the  assembly  and  retention 

of  materials  into  a ROVAC  vane  structure. 

2.  Through  the  use  of  the  advanced  vane  materials  and  structural 
design  changes  analyzed  in  this  program,  total  weights  and  operational 
bearing  loads  were  reduced  by  3S%.  In  addition,  the  use  of  advanced 
composites  allows  for  a design  optimization  toward  the  thinnest  structure 
possible  from  either  a manufacturing  or  operational  point  of  view. 
Specifically  the  analysis  considered  vane  sizes  down  to  40%  of  original 
size.  For  the  reduced  thickness  cases,  the  vane  was  structurally  sound 
(in  the  analysis  plane)  and  total  weights  and  operational  bearing  loads 
were  reduced  by  66%.  Both  the  35%  and  66%  reductions  can  also  be 
accomplished  while  concurrently  reducing  the  peak  vane  stresses  by  30% 

to  60%. 

3.  The  primary  objectives  of  the  study  were  to  reduce  the  amount  of 
tip  and  edge  deflections  occurring  within  the  vane  structure.  Based  on 
the  results  of  this  study  it  was  concluded  that  deflections  can  be 
significantly  improved  through  materials  and/or  design  variations.  Within 
the  analysis,  the  most  significant  changes  occurred  from  the  composite 
materials  substitution.  Out  of  the  nine  materials  variations  considered 
within  the  analysis,  the  two  best  showed  tip-deflection  reductions  of 

68%  and  59%,  and  concurrent  edge-deflection  reductions  of  72%  and  70%, 
respectively.  The  removal  of  vane  axles  and  retention  boots  increased 
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the  improved  deflections  by  approximately  20%.  These  increases  are  only 
minor,  however,  since  the  original  materials  changes  would  still  result 
in  62%  and  50%  improvements.  The  effects  of  reducing  the  vane  thickness 
were  generally  positive,  in  that  due  to  the  decrease  in  inertia  effects, 
deflections  actually  decreased  even  though  the  structure  was  thinned 
down  substantially.  Actual  values  were  only  minor,  however,  with  a 
maximum  9%  improvement  being  noted  in  the  0.4t  case. 

2.  RECOMMENDATIONS 

Based  on  the  results  of  this  program,  it  becomes  evident  that  further 
study,  analysis,  and  testing  is  required  to  fully  optimize  an  advanced 
vane  configuration.  In  this  regard,  the  following  considerations  or 
recommendations  are  offered: 

a.  Materials  Consideration 

Based  on  initial  ground  rules,  the  analysis  discussed  in  this  report 
was  restricted  to  advanced  graphite  composites.  These  types  of  materials 
have  several  basic  disadvantages,  consisting  of  high  friction,  abrasiveness, 
cost,  and  potential  corrosion  inducements.  In  addition,  from  the  analysis, 
it  was  noted  that  the  primary  advantages  of  strength  were  not  being  fully 
exploited.  By  sacrificing  strength  and  perhaps  some  stiffness,  the 
associated  materials  disadvantages  might  be  offset  through  alternate 
materials  choices.  Alternate  materials  choices  offering  improvements 
in  these  areas  include  Nomex  and  Kevlar,  both  of  which  could  readily  be 
incorporated  into  a ROVAC  vane  structure  as  either  a compression-molded 
or  laid-up  structure.  Both  materials  would  eliminate  corrosion  and 
cost  problems  at  the  price  of  reduced  but,  most  probably,  adequate 
strength.  In  addition,  these  materials  offer  further  improvements  in 
the  area  of  weight  and  friction.  By  further  reducing  the  strength 
requirement,  substantial  improvements  could  be  obtained  (relative  to  cost 
and  friction)  through  the  use  of  fiberglass  and/or  powdered  additives  to 
the  base  resin  system  of  the  vane  structure.  This  friction  improvement 
may  contribute  to  the  reduction  of  internal  heat  generation  within  the 
ROVAC  machine. 
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Based  upon  the  previous  discussion,  four  very  important  recommendations 
can  be  offered. 

1.  A preliminary  structural  analysis  should  be  conducted  defining 
the  penalties  and  payoff  for  various  alternate  materials.  Materials  to 
be  considered  should  include  Nomex,  Kevlar,  fiberglass,  and  woven  Teflon 
fabrics.  In  addition,  resin  system  additives  and  supplements  should  also 
be  considered  within  the  analysis. 

2.  It  is  recommended  that  the  analysis  discussed  above  be  conducted 
by  simply  varying  the  MAT2  properties  of  the  NASTRAN  program  contained  in 
Appendix  C of  this  report. 

3.  For  the  recommended  analysis,  physical  properties  of  some 
recommended  materials  may  have  to  be  generated.  It  is  recommended  that  a 
materials  test  and  evaluation  program  be  initiated  to  define  friction, 
wear,  strength,  and  deflection  properties  for  potential  ROVAC  vane  materials 
for  which  no  current  data  exists. 

4.  Many  similar  or  even  identical  materials  may  be  considered  in  a 
current  FEM  advanced  composite  bearing  effort.  It  is  recommended  that 
coordination  be  maintained,  with  these  efforts  as  they  may  prove 
invaluable  in  the  selection  of  an  advanced  composite  ROVAC  vane  material. 

b.  Analysis  Verification 

Although  the  reported  analysis  agrees  with  ROVAC  deflection  and  loads 
data,  it  is  recommended  that  a correlation  test  program  be  accomplished 
which  would  demonstrate  the  complete  validity  of  the  AVM  and  NASTRAN  models 
generated  in  this  report.  The  correlation  program  would  not  have  to  be 
elaborate  and  could  simply  consist  of  a representative  static  loading  on 
a baseline  ROVAC  vane  noting'overal 1 deflections  and  several  selected 
strain  values. 
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c.  Design  Variations 

All  of  the  designs  considered  within  this  study  were  restricted  to 
the  existing  ROVAC  design  (Figure  4).  For  this  design,  geometric 
variables  were  constrained  to  values  which  may  not  necessarily  be  an 
optimum  from  a materials/structural  relationship  point  of  view.  It  is 
therefore  recommended  that  geometric  variations  from  Figure  4 be 
analyzed  for  more  optimized  materials  selection.  These  optimums  could 
then  be  compared  to  optimums  from  an  air-cycle  machine  point  of  view  and 
trade  studies  initiated  for  a fully  optimum  approach.  These  design  runs 
could  be  very  easily  accomplished  through  the  use  of  the  AVM  program 
developed  in  this  report.  In  addition  to  the  geometric  changes,  the  AVM 
program  also  can  rotate  the  materials  axis  to  any  value  in  four  distinct 
areas.  Applications  of  this  capability  may  lead  to  still  further  design 
improvements. 

d.  Refined  Analysis 

All  of  the  design  changes  studied  in  Section  V assumed  that  some 
form  of  bearing  support  pin  could  be  structurally  integrated  into  the 
advanced  composite  region.  Currently,  several  approaches  exist  for 
accomplishing  this  but  all  would  require  modification  in  areas  of 
maximum  stress.  In  view  of  this  fact,  it  is  recommended  that  these  areas 
be  refined,  and  the  actual  attachment  method  be  modeled  to  determine 
the  best  method  for  bearing-load  transition  into  the  vane  structure. 

Here  again,  the  AVM  program  could  be  applied  so  that  it  possesses  the 
capability  to  readily  refine  any  area  within  the  vane  structure. 

e.  Margin  of  Safety  (MOS)  Study  (Recommendation) 

Although  the  stress  levels  noted  in  this  report  appear  low,  an  MOS 
analysis  should  be  conducted  to  pinpoint  relative  vane  strength.  This 
analysis  would  then  enable  realistic  projections  to  be  made  relative  to 
the  degree  of  strength  reductions  permissible  for  the  advanced  materials 
previously  recommended. 


81 


AFFDL-TR-77-49 


appendix  a 


PROGRAM  VANE  LISTING  IN  AVM  FORTRAN  CODE 


I 
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APPENDIX  A 

PROGRAM  VANE  LISTING  IN  AVM  FORTRAN  CODE 


i??SCEDI\0  Page  3L.aNK_N0T  FILMED 
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1 


«7XK.S»StA.U'(»IlX-.  07tai5i^SCHiUUXi»«vl4. 

FTM. 

LGO. 

PFWIND,TAPt7, 

COPySRF. TflPE7, OUTPUT, 

PfWI Nn,TAP£7 . 

_ £aPTilF^TAP£Z»PUHCH.-  . . _ _ 

EXIT,';  . 

PP.OGRAH  VANE  ( I NPUT  , T APE5=  I NPUT , OUTOUT,  T APE5=  OJ TP UT  , T APE7) 

r 

C THIS  PROGPAH  P.FADS  X ANO  i DATA  FROH  3AROS  AN3  CONSTRUCTS 
C A GRID  FOP  OUTPUT  IN  NASTRAN  PORMAT 

rOHM09  /VANX7X<700) ,X(700» ,XP(5),XTH{6»,NO,TH<V) 
C0HM0N70UM3AR7XR001,yR001, YR002,NR00,NG1,MG2,IT3T ,NR1,NR2 
C 

C RIAD  CARD  INPUT 


CALL  RIACIN 

_ C - 


c 

CALCULATF  X-Y  GRIO  ANO  print  GRID  CAROS 

CALL  GRID 

r 

c 

CALCULATE  ANO  PRINT  VALUES  FOR  CQU432  CARDS 

C.  - 

CALL  COUAOP 

c 

CALCULATE  CPAR  VALUES 

IF  (NROD.NE.O)  CALL  COAR 

c 

. c 

CALCULATF.  SINGLE  “OINI  CONSTRAINTS 
call  SPC 

r, 

f 

CALCULATP  HATFRIAL  ANO  PROPERTY  DATA 

EALL  SALEIN 

ENOFILE  7 

PTCP 

END  . . . 

c 

r 

SURROUTirE  REAOIN 

PEAOIN  READS  TRP  pOLLORINE,  VALUES-- 

r 

c 

1,  roL  1-5  - npx  - no.  of  x va^jes 

CARD  1,  COL  &-iQ-  MPY  - NO,  OF  Y VALUES 

c 

c 

p,IFMT,ANP  XJ  A'’E  REAO  in  SPECIAL  FORiAT 

ALL  OTHER  VALUES  READ  HITH  6F1 0. 0 FORMAT 

COMMON  /VANX/X  (700  ) ,Y(700),X’(5),XTH(5>,N?,rH(:*) 
rOMMON/OUMOAR/XROOl ,YR001 , YPOO2,NR0n,NGl,NG2,ir0T 
CQMMON/AGAIN/XlNl7b» ,YlNt7P) , NPX ,NPY,NEL 

,NR1,NR2 

COMMON/E XT R A/RAD  IDS ( A > , TH I CK ( 4) , X J ( A ) , E ( 8) , XM U ( 3 ) 
1 ALDAIAtS) ,XI (2,4) ,VEL, AXIS.GRAV 

pIMpnSIO)  TR(2) 

LOGICAL  PRINT 

,RH0(8)  , 

EFAO  I5,ED3)  NPX, NPY, PRINT 

£,30  pOphaT  (2I5,L1) 
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asAD  ixiNui , 

fitan  ,5,Fio>  (Yi-Mi) ,i=i,Npy) 

51C  rOR'IAT  (8Fia.  0» 

PFftn  (F,Fi;)  yo 
RtAO  (SftlOl  XTH 
PFAn  (5,  MO)  TM 

KEAa  t5»t201  .'«0Q»XR0Dl,YiOCl^XR002»)il£O?II« 

F?C  FORMATdn.YFlO. 0,110) 

YS(l):YRODl 
YR(2)  = YRrn'> 

K’PYM1  = MPY-1 

IF  (NEOPIIN.GE. 1)  GO  TO  30 

IF  XAlROn.EO.a)  GO-  TO  3C  

DO  15  I=l,NROn 
DO  10  .I=1,M°YM1 

IF  (YIM(J)  .LT.YR(I)  .ANO.YINI  )*l).3T.fR(I))  30  TO  5 
GO  TO  10 
5 HPr  = H»YH 

JN-MPX  . .. 

7 YIN (JN) =YIN( JN-1) 

JN=  JN-1 

IF  (JN.GT.JH)  GO  TO  7 
YIN  (Jf  1)  :YR(I) 

IP  rONTINUF 

IE  CONTINUE  - 


RFAOS  IN  All  of  TME  materials  and  PROPERTIF-;  ffPF  qaTA 
30  CONTINUF 

PEAO  (5,100)  IF,)T 

IF  tlFMT.NE.U  GO  10  ILO  . - 

PFAO  (5,11,0)  XJ 
SAC  FORMAT  (l,CX,',F10, 0) 

PEao  (5, MO)  THICK, RADIUS 
RFAO  (5,f30)  E 
53C  F0R“AT {8(?X,A8) ) 

READ  (S,Fici  xMU,RHO,«EL,AXIS,GRA)/.  - . . 

pfaO  (5, Fin)  ALOATA 
IF  (PRINT)  CALL  LISIIN 
RETURN 

IOC  PRINT  600 

6CC  F0R“aT(*  THIS  AR^A  NOT  I M® LE I ENTED . • ) 

STOP 

FMO 

SUBROUriNF  GRID 

THIS  ROUTINE  calculates  THE  X-Y  GRID  POINTS  f^qn  THE  INPUT  l/ALUES 

CQHHON  /V4NX/X  (700)  ,Y(7001  ,X-®  tSl  ,XIH(5i,NP^lHUl  

COMM0N/ArAIN/XIN(7E) ,YIN(75) ,NPX,NPY,NEL 

COMHON/OUMgAR/XRODl ,YR0D1, YROOZ,NROn,NGl,NG2,IT3T ,NR1,NR2 
NO=  0 

on  15  j=i,npy 
oo  15  K=1,NPX 
NPiMPH 
X(NP)fXI)'(K) 

Y(NP) sYlN( J) 

15  CONTINUE 
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IF  (NROO.tO.o)  r,o  TO  5CC 
OALOULATF  '■.oIO  valufs  for  fir^t  bar 

XINC-ABS (XIMCl) -KROOl) /5. 

. . - 

PO  lOJ  J=l,5 
NGl^NGlU 

X<NG1)  =XFO01*-(.)-l)*XINC 
YlNGll =VROni 
irc  rONTI'J'JE 

00  I MAVF  TO  CALClJLATt  NFW  POINTS  THROUGH  5TRJCTJRF? 
NR1  = C 

no  110  J=1,N‘>Y 

IF  (YRODl.EQ.YIN(J)  ) GO  TO  HO 
lie  rONTIHUF 

MO  MAIGH  ..  - - - 

HAY'"  NFW  FRIO  POINTS 
no  120  J=1,NPX 
Nr.l  = NGlH 
y(NGl) =XIN{J) 

Y (NGl) =YF001 
12e  CONTINUE 
GO  TO  IBC 
lie  t'Rl  = J 

Ron  2 
MR  2 = 0 

i£C  IF  tNaQO.LT.21  QQ  TO  SC-0  - 

no  170  J=1,NPY 

IP  {YROnF.rO.YlNtJ) ) GO  TO  230 
170  PONTINUE 

NO  MATCH  - NEW  GRID  POINTS 
NG2=N31 

no  2'0  J=1,NPX 
NG2=NG2h1 
X(NG2»  =XIN  ( J) 

Y(NG2I =YPnD2 
70C  PONTINUF 
GO  TO  SOI 
250  NR7=J 

HFI7F  grip  cards 
50L  ITOTriAxr (NO.NGl ,NG2) 

WRITC  (7,7331 (J,XtJl,YtJ),J=l, IT OTl 
70  0 FnRMAT(fcHGRn,4X,I8,HX,2F8.2,9H  0.00) 

HEIUSN 
I Nn 

SUBROUTIie  COUAD? 


calculate  AND  PRINT  COUAD?  CARDS 
POMMON  7VANX/X  t730),Y(700),X7(5),XTH(5),N’,TH(Rl 
common/ AGAIN/XIN1751 , YIN ( 7 E) , NPX.NP Y, NFL 
PIMFN3I01  IGP(4) 

NEL=C 

NXL=NPy-l 
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nn  50  JsltNYL 

IGP(ll=IJ-l)*NPX  _ _ 

no  30  K=1,NXL 

N£L=M£LH  

lGP(i»=ir-p(i»n 

1GEL121  ilGP  I U-tl 

IGP<3)=ir.P«?)fNPX 

IGPfl.)  =IGP(1J«-NPX  

FIND  PID 
no  5 M = l,c» 

IF  -IX^tXGF  t211^at*XP^rtJ.4N3.XII&PI2iJ«LE.X5i'tHXl-j;0  ICL-7 

5 CONTINUE 

*»  SET  I‘>=0  SO  E990R  CAN  BE  CHECKED  (THIS  SHOULD  NEVER  HAPPEN) 
IP=0 

GO  TO  10 
7 IP=H 


FIND  TH 

10  PO  15  H=l,4 

IF  <X(I6P(7)1  .GT.XTH(H) .ANO.X  (IGP(?)) ,LE.XTH(Hfl) ) GO  TO  2D 
15  CONTINUE 
•*  00  SANE  AS  ABOVE 

IHX=£. . - . - - 

GO  TO  25 
20  THX=THtM) 

25  WRITE  (7,700)  N-^L,  IP , IGP, T MX 
70C  FORMAT (6HOOUA02,2X,6I8,F5. 3) 


30  CONTINUE 
5C  CONTINUE 
RETURN 

IMP - . 

FUBROUTIf'F  LI'^TIN 

LIST  WRITES  ALL  INPUT  PARAMET' RS 

COMMON  /VANX/X ( 7 00 ) ,t ( 7 C 0 ) , X” (5) ,XT HI 5) , N’ , TH ( 4 ) 

XaM.>lDN/OU-MaAR/Xfi001«YRODl»TgD-D2«Jiiafla»AU;i»NS2,II3T»NRl,NR2  

C0MM0N74GAIN/XIN(7E)  ,iriN(75)  ,NPX,NPY,NFL 
C0MH0N/CXT°A/RADIUSI4) ,THICK(4) ,XJ(4) ,E (8)  ,XMU(8)  ,RHO(e)  , 

1 AL0ATA(8) ,XI (2,4) , VEL, AXIC.GRAV 
WRITE  (6,600)  NPX, IXIN(I) ,I=1,NPX) 

60C  FORMAT (’IPPOGRAM  re  AO  *,I3,»  X-VALUcS  AS  FOLLOWS--*/ 

1 (lOX, 10F19.2) ) 

WRITT  16,610)  Npr, (riN(i) , i=i,NPr) 

610  FORMAT  (///5X,»  AND  READ  *,13,*  Y-VALUES  AS  FOLLOWS—*/ 

1 (lOX, 10F10.2)  ) 

WRITE  (6,620)  XP 

62C  F0RMAT(///5X,*PROP£RTY  ID  OIV ISIONS--* ,5Fie . 2) 

--  WRITE  (6,630)  XTH,TH - 

63C  FORMAT (///FX, 'MATERIAL  PROPERTY  ORIENTATION  ANOtE  INFORMATION — •/ 

1 10X,*OATA  DIVISIONS--*, 5F10. 2/, i:X, 

2 'ORIENTATION  ANGLES--*, 4F1 0. 2) 
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WRITF  (6,b<*a)  XR0Dl,YR001,YRan? 

bUC  FORMAT (///5X»»XR001  = *tF8. 2,*,  YROU 1=» , F«. 2 , * , YR002=*,F».  2 > 

HRITF(6,f 991 

699  FORMAT  (///FX,»PROPFRTY  1 2 t - -4  - 

X F 6 7 8 * » 

URTTF  fft-Ynni  y.i,TnTr.XTRAnTii>; 


70 C FORMAT (/7/5X,»J=*,42X,4F10.3//5X,*THICK  = *,l*F10,  T// 

1 5X,*F,ADIUS=f,37X,4F10.3»  __ 

HRITi^  (6,7101  F,XMU,RHO 

71 C FORMAT (//5X ,*E=* ,4X , SA1C//5 X, *XMU=* ,8F10.3/75X, *RHO=*, 8F10, 31 
WRITE  (6,7201  VFL, AXIS, GRAY, ALOATA 


1 *,  rRAV=»,F10.3//5X,*ALOATA=*,8F10.3) 

RETURN 

ENn 

SUBROUTINE  COAR 

C 

X CAinillATFS  AMO  3RTNTF  VAIIIFR  FOP  OOAR  CAROS 

SINCERE  APOLOGIES  TO  STRUCTURED  PROGRAMMING 


C0MMnN/AnAIN/XIN(76» , YIN ( 75) , NPX ,NP Y, NEL 

COMMON  /VANX/X (700) ,Y(700)  ,XP(5) ,XTH(5) ,NP,TH(4) 

rnMM0N/0UM9AR7XR0Ol,YR0ni,  YROD2,NR3n,NGl,Nr,2,  IT9T,NR1,NR2 

OIMENSTnK..IGP(7) 

NFX(NN)=  (NN-l)*NOXH 

r 

c ralculatr  all  horizontal  cpar  values  for  ROO  1 

IPIO=5 

X1=0. 

X2.?l. 

X3  = 0. 

IF=1  


N^LO=NEL 
IGP (1) =NF 
no  50  J=l,4 

. --tlELa^NELF-fcl 

iGP(i)=ir.p(i)»i 

IGP(2)=IGP(l)fl 

WRITF  (7,700)  NELB,IPI0,IGP,X1,X2,X3,IF 
7C0  FORMAT (AHCOAR, 4X ,41 8, 3F8. 2 , 1 9) 

EC  rONTINUF 


WHERE  DO  I GO  FROM  HERE? 
iPin=7 

IF  (NRl.PE.O)  GO  TO  100 
no  70  J=1,NPX 
IGP  (1)  =IGP(1)H 


r 

X.- 


NELB=NELRH 

WRITE  (7,700)  NELB,IPI0,IGP,X1,X2,X3,IF 
70  CONTINUE 
GO  TO  12C 

-FIHO  OLD  Gftia  POINTS 


100  IGP(1)=IGP(1)H 

IGP(2)  =(NR1-1)*NPXH 
NFLO=NFLP«-l 
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WBITC  <7,7-004  NELBtIPIO,IgP,K-ttX-2»X^rIg 

I6P(H=irp(a)-i 

_ NPXMlsNPX-1— 

no  110  J=1,NPXH1 

IGPCH=ICPll>H  

ir,P  (2t=IGo<lJ»l 

WfcLBsNfLRKl 

MRITF.  (7,7C0»  McLB,IPIO,ISP|Xl,X2,X3,IF 
_ 110  CONTINUE  _ — _ 

r 

C _ 00  I H*WE  A SECONO  ROO?  

120  IF  (NROO.LT.2>  1,0  TO  200 

I«»ID=B 

IF  (MR2.CO.01  GO  TO  130 

IGPIll=tNR2-l)*NPX  

GO  TO  14C 

_ 110  IGPIllsNEl  - 

140  NPXM1=MPX-1 

no  iso-j=i.»MexMi 

IGPdl  =IPP(H»1 

_ IGP(2l5lGPClJH - 

f4ELB=MELP«-l 

WRITE  C7, 7001  MELB, IPIO,IGP,Xl,X2,X3,IF  - - 

150  CONTIXUe 

C 

r NOW  FOR  THE  VERTICAL  BARS  ON  ROD  1 

_ C find  tie  POINT  

20C  Xt=l. 

X2r0,  _ 

IPI0=6 

IF  tMRI.EIUJt-  GO -TO  250 

NUPsMINO  <3>MPX-HR1» 

NDNsMINO  l3fNRH  

C 

C CRAR  CAROS  FOR  EXISTING  GRID  POINTS  

NF.LB^NELO*! 

IGP  Cl>  s-< KR-1--2J  » N°X » 1 

IGP{2»=NP»5 

HRITE  (7,700*  NELB,IPI0,IGP,X1,X2,X3,IF 
IGP(1»=IGP(2» 

_ IGP(2»=NPl*NPXf 1 

NEL8=NELP«-1 

HRIIE.  (7,7-001  NELB,IPI0,IGP,X1,X2,X3»1F 

C 

NT0sNUP*N0N  

IGP(1>  s(NRl-NON-l»*N»X*l 

— DO  210  J=l,NTO  

NELB=NeLP*l 

lGP(2>aIGP(m-NPX 

HRITE  (7,700*  NEU8,IPI0,IGP,X1,X2,X3,IF 

IGP(l) aIGP(2)  

210  CONTINUE 

RETURN  

C 

— C C-BAR-IUfiOS  FOR  NEW  GRID  POINTS 

250  00  260  J=2,NPr 

IF  (TIN(J)  .GT.rROOl.ANO.riNtJdl.LT.rROOl)  SO  TO  265 
260  CONTINUE 
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?&.£.  MlOPaJ 

NROT=J-l 

NUP=MIN0  (2*M>>ir-NT0P>  

NnN=MIN0  <2,  CNPOT-in 

IGP (1» =MFX(N30I-N0N)  

27C  IGP(2)  =ir,P(l)fNPy 

WlR=Mt:LR»1 

WRITE  (7, 7001  NPLB,I“IO,ISP,yi,X2,X3,IF 

IGP (II =IGP(2)  

IF  IIGP( 1) .EO.NFXINBOT) I GO  TO  280 

GO  TO  27C  

2*0  ir,P(2)=NPF6 

HFLRr«|rLPt.1 

WRIT'^  (7,700)  NFLB,  IPIO,IGP,Xl,X2,X3,IF 

IGP(1I=IG“(2)  

IGP(2)=NFX(NTOP) 

NELfi=NELFH  

WRITE  (7,700)  NELB,IPID,IGP,X1,X2,<3,IF 

290  IGP.I1I=TGP(?) 

IGP(2)=IGP(1)f*(PX 
NELB=MELP«-1  

WRITE  (7,700)  NELB,I'>I0,IGP,X1,X2,X3,IF 
IF  (IGP(2).LT.NFX(NUP4'NT0P))  go  TO  290 
C 

C OO  TWO  STFANGF  OIAGOMAL  BA.R5 

I'»IO=5 

NELB=yFLPfl 

WRITE  (7,700)  NFLB,  IPIO,NFX(tJOOT)  ,'IPF‘»,X1,X2,X3,  IF 
NELB=‘)£LPH 

WRITF  (7,750)  WELB,IOIO,NP*5,NFX(WTOP),X1,X2,X3,IF 


RETURN 

END 


FUBROUTIF'^  CALOIN 
C 

C this  ROUTINE  CSLC'JLATF.S  HOTERIftLS  AND  PROPERTIES  TYPE  DATA 

_ C 

rONHON  /VANX/Y  (700)  ,Y(700)  ,XO  (5)  ,XTH(5)  ,N»,rH(!,) 
CO*1NON/AGAIM/XI)4(76)  , YIN ( 75)  , NPX,NP Y,NEL  - 

COMMOM/DUH«AR/XR001,YR001, YRO02,NROD,NGl,NG2,IT0T ,NR1,NR2 
COHNON/EXTRA/RODIUS (4) , THICK (4) , XJ ( 41 , E ( 0)  , XNU ( S) ,RH0(8) , 
1 ALOATA(S),XI(2,4) , VEL , A XI S, GRAY 

niHFNSIO).'  -A-C4) 

PAT  A »I/3.  14159265/ ,7ERO, ONE, ON£NE570.,1. ,-l.  / 

C 


no  1C  .1=1,4 
A(J)=PI*f A0IUS(J)**2 

XI (l,J)=XI(2,J)=oi*RADIUS(J)**4/4. 

nONTIAIUE 

C 

WELC=tfEL/6Q. 

IF  (GRAV.EO.O)  GO  TO  699 

GACF.L=GRAV 

GO  TO  701 

699  CACEt=ALDATAtH»ALOATA-t2T 

701  CONTINUE 
C 

I r WRITE  EXTRA  CAROS 
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70C  FORMAT <*PQUft02  *t21fl>F8.4> 


no  115  J=l,4 

JSsJtV  - 

V91TF  (7,705)  JS , JS , A( J) ,XI ( 1 , J) , J, J, XI ( 7, Jl , X J ( J) 

-ZJ5-  F-OBtUT  (aHPOAR* ,2I16,2E16.7,»QABeP£F-»,ll/7HTA9C01 

115  CONTIIUF 

_ DO  120  J=1,S . 

HRITf  (7,710)  J,5(J) ,XHU(J) ,RHO(J) 

71C  FOBMAI  {•(‘ATI  _ *,18,  A A , SX,  2F  0,4)  ^ 

120  contimhf 


IOS=NPX*M(»T*-l 

IT0TP1  = IT0TH  _ 

HRITF  (7,715)  IRS,IT0TP1 

715  FORMAT  (*5:PC1*,9X, *500  123466*,2IA)  

HRITF  (7,717)  ITOT 

717  FORMAT  f *SPC1*,-9X-,»501*  , 7X  , »6»  ,7J(«t.l THRUF,I») 

IBSlsIRS^l 

IBS2=iaS  + 2 

HRITF  (7,720)  IBSl,Ii3S2 

72C  FORMAT  (•?.PC1*,9X, *602*, 7X,»2* ,218)  


725  FORMAT(7PR'^ORCC*,15X,*20*,I15,l6X,E16.7,*aRFORSEl*/ 

1 aH»RF0RC£l,  3F8.4) 

HRITE  (7,7T0)  GACEL , ZE RO , ONEMEG , ZERO 
730  FORMAT (5HGRAV* , 17X, * 21* , 16 X,E 16. 7, F1&. 4 , ’aGRAWl*/ 

1 6H)’GFAV1,2X,  2F8.4) 

MRITE  17,7351  -OM£,OMc,ONE 

735  FORMAT (»PAPAM  HTMASS  . 0 0 2» e8*7*“ARAM  GROPNT* ,9X,*0*/ 

1 *LOAD*,10X,*30*,2Fa,4,6X»*20*,F8,4»5X,*21*)  

YAX=YIM(1)«-AXIS 

HRITE  (7,740)  ITOTPl,ZtRO, YAX,7ER0 
74 C FORMAT(*rRIO*,4X,ie,8X,3F8.3) 

HRITE  (-7, 741) 

741  FORMAT (•INOOATA*) 

RETURN  _ 

FNO 

SURROUTIKE  SPC 


C0MM0N/AGAIH/XIN(75I ,YIMt75) ,NPX,MPT,MEL  

COMMON  /VANX7X(700) ,Y(700) , X® (5 ) , XT H(5) ,N»,TM(4) 

C0HH0N/OUMBAR/XROI)l,rROOl,yR3O2,NR009NGl,NG2,ir0T,NRl,NR2  

OIMENSIOl!  ISRlSu) 

-ISP(l)a99 

HL=0 

ICON=15  - 

N$PC=NPY 

DO  10  J=1,NPY  

ISP(JH)  = ISP(J)  n 


HRITE  (7,700)  I9P(JY1) ,IC0N,NL 
700  FORMAT (4HSPC1,4X, 318) 

1C  f'ONTINUF 
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IF  IMRl.NE.O)  r,0  TO  15 
NSPCsNSPCH 

ISP(NS»CH)=ISP(NSPCIH 
HRITE  (7»7Cat  ISP  tNSPC*ll , rC0M,N61 
15  IP  <MR?.T'E.  0>  00  TO  20 

NSPCaMSPr.H 

ISP(NSPC«-1»=ISP«NSPCI  fl 

MPITE  (7,7001  ISP(NSPC*11 ,IC0N,NG2 

WRITE  SPC  TOTAL  CARO 
20  NSPC=NSPr«-4 

MT.ARDsHSPR/a 

IF  (MOnCNSPC,8l  .NE.  0»  NCARDr'^CAROU 

ISP(HSPC-21=5DB 

ISP(NSPC-1)=501 

ISP(NS?C)sS02 

NCONT=10 

IF  (MCARP,EO»l)  GO  TO  20J 

WRITE  (7,710)  (ISP(I),I=l,8),NCONT 
71C  FORMAT (6HSPCAD0,2X,8IB,3HA3C,I2j 
IF  (NC1RD.LF.2)  GO  TO  100 
NCR=MCARC-1 
11=1 


Il=Il«-8 

IS=Il«-7 

IF  (IS.GF.MSPC)  GO  TO  70 

IO=NCONT 

MCONTsMCONT*! 

WRITE-.  (7,7201  IQ,  ( TSPtll  , T si  1 ,IS: 

50  CONTINUF 

ii=ii»a  

70  IS=NSPn 
I0=NC3NT 

720  FORMAT(3M»BC,I2,-^X,BI6,3HARC,I2) 

73C  FORMAT  (3H»-0C, 12, 3X,  818) 

PETURW 

100  11=9 

IS=NSPC 

IQg  10. 

NCONT=ll 
GO  TO  75 

200  WRITE  (7,71.0)  (ISP(I)  ,I  = 1,NSPC1 
740  FORMAT (6HSPCADn,2X, 818) 
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appendix  B 


TEST  CASE  PLOTS  AND  ABBREVIATED  LISTINGS 
FOR  TEST  CASES  1 THROUGH  4 
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0 . uC 
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COUAO^ 

3 

2 

3 

. 

18 

17 

0.  JOJ 

CQUA02 

“♦ 

2 

5- 

19 

l3 

0.  0 00 

CQUA02 

5 

3 

9 

20 

19 

0. 3C0 

C0JA02 

0 

3 

0 

7 

21 

2u 

J.  0 30 

CQUAD2 

? 

7 

0 

22 

21 

0.000 

CQUA02 

8 

& 

9 

2 2 

22 

0.000 

C0UAb2 

9 

4 

Q 

10 

2s 

23 

0.000 

C0UA02 

10 

t 

10 

11 

2-j 

2i. 

u . 0 G 0 

CU0AD2 

11 

11 

12 

2e 

25 

0 . 3 CO 

CUUA02 

110 

s 

124 

l20 

139 

lio 

0 . 9 uO 

CQUAD2 

117 

s 

126 

12c 

140 

139 

0. 0 90 

CBAR 

118 

J 

141 

142 

u . 0 0 

l.GO 

C.OO 

1 

CBAR 

119 

0 

l‘-2 

ls3 

0 . L u 

1.90 

0 . b 0 

1 

CBAR 

120 

5 

ls3 

14S 

u.  CO 

V • L ^ 

0.00 

1 

CBAR 

121 

5 

144 

I4i 

0 . C u 

1.90 

w . U J 

1 

CBAR 

122 

7 

IsS 

16 

G . 0 0 

1 • J ij 

G..00 

1 

CBAR 

123 

7 

16 

lb 

0 . C 0 

i.  20 

0.00 

1 

CBAR 
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7 

1C 

17 

b . 0 0 
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0 . 0 
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17 

13 
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7 

18 

19 
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1 

CBAR 
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7 

19 

20 

b . 0 0 

l.OJ 
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1 

CBAR 

123 

7 

21 

b . b 0 
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1 

CBAR 
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7 

21 

22 

0.00 
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O.CC 

1 
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130 

7 

22 

23 

b.  GO 
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L • J u 

1 

CBAR 

131 

7 

23 

2s 

O.CO 
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0 . 0 J 

1 

CBAR 
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7 

2s 

25 
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1 

CBAR 

133 

7 

26 

CO 

C.  00 

1 • 0 u 

e . 00 

1 

CBAR 

134 

7 

20 

27 
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1 • o 0 

Cl  « J d 

1 

CBAR 

13? 

i 

27 

28 

0 . u b 

X • w 0 

C « d w 

1 

CBAR 

13b 

3 

o6 

8o 

0 . b 0 

1.  jO 

O.-Ou 

1 

CBAR 

137 

3 

oo 

87 

G . C 0 
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1 

CBAR 

138 

8 

67 

08 

b . 1 0 
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u • 00 

1 

CBAR 

139 

J 

to 

So 

b . C b 

1 • w 0 

0 • Co 

1 

CBAR 

IhO 

0 

89 

9i, 

O.CO 

1,'30 

l.OO 

1 

CBAR 

l4l 

8 

90 

91 

0.00 
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b.  UQ 

1 

CBAR 
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0 
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92 

0 . u 0 
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8 

92 

93 
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3 

94 

95 
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1 
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C84R 

l-*b  t> 

95 

9b 

0.  oC 

1.00 

u.OD 

1 

C8AK 

i‘t7  i 

9d 

97 

0.00 

l.uO 

Q.  OJ 

1 

GEAR 

1h4  a 

97 

95 

0.00 

1 . w 0 

D « Ou 

1 

CBAR 

l-^B  j 

1 

1-3 

1.00 

0.00 

C • C 3 

1 

GBAK 

15J  0 

Ist 

23 

l.CO 

0.00 

b » uO 

1 

CBAR 

131  3 

-13 

1 

1.00 

0 . u 0 

0*GJ 

1 

CBAR 

152  3 

1 

l3 

1.00 

0.00 

0 • u u 

1 

CBAR 

153  j 

15 

29 

l.CO 

0.03 

0 • uC 

1 

CBAR 
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29 

*.3 

l.CO 

3.03 

0*00 

1 

GBAK 

153  0 

h3 

57 
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0.00 
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1 
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It 
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28 
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^2 
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7 0 
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Figure  B-5.  Test  Case  2 (Elements) 


Figure  B-6.  Test  Case  2 (Nodes  and  Elements) 
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Figure  B-7.  Test  Case  3 (Nodes) 
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Figure  B-8.  Test  Case  3 (Elements) 
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Figure  B-9.  Test  Case  3 (Nodes  and  Elements) 
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Figure  B-10.  Test  Case  4 (Nodes  and  Elements) 
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VANE  DEFLECTED  SHAPES  FOR  22  ANALYSIS  CASES 
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